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Chapter 1 


INTRODUCTION 


Rotorcralt noise has become a subject of considerable inten'st to the aerospace com- 
munity. Manulact urers ol civil rotorcralt arc facing increasingly stringent federal and local 
noise* regulations, and military aircraft makers are* being required to reduce noise* signature's 
in the* interests of battlefield survivability. 

A compute*r program called ANOPP has been developed! at the* NASA Langley Research 
(enter for prediction of the noise* r**coi ve**l at the ground due to aircraft. The* ANOPP 
program has evolve**! into a system with threx* subsystems: a system for CTOL aircraft, 
the Propeller Analysis System (PAS), anel ROTO NET. The* ROTO NET subsystem was 
developed to predict the external noise generated by a given lielicopler flying in a specified 
flight condition. Ibis program calculates several noise' sources on the helicopter, predicts 
the propagation ol the noise to farfield observers, sums the various noise components, and 
computes various noise' metrics for use' in certification studies. 

This elocument pre>vieles intmeluctory information fe>r newice users e>f th<‘ ROTONKT 
code'. Hrie'f re'vie'ws are' preseuiteel on; acoustics thee>r\\ rote>rera(t ne>ise' phenomena, ro- 
tonkt cejncept.s, input data cre'at ie>n. <x)cU' e*xec utie>n. anel interpretat ie>n of results. In ad- 
dition, ROTONKT users will alse> need the fbllenving documents; “Aircraft Noise' Ptveliction 
Pre)gram Theoivt ieal Manual [1], “Aircraft Neiise Predictie)ii Program Pser’s Manual" [2], 
“'The' NASA Aircraft Noise Preelict ie>n Pre>gram Imprenvel Pmpelle'r Analysis System” [2]. 
anel "Aitrraft Ne>ise‘ Preelict ion Program J 'hoe)rot ical Manual" [I]. An aelditional reference 
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which can prove helpful for developing non-level flight trajectory predictions is “An Intro- 
duction to High Speed Aircraft Noise Prediction" [5]. Those documents provide detailed and 
comprehensive information on the theoret ical methods used, t, he input and output quantities, 
and the ANOPP Control Language' syntax. 

Since theoretical and programmatic quantities are discussed in this work, several type 
fonts are used to distinguish between forms. All ANOPP input and output quantities are 
printed in Courier font . for example t he SFIELD(PHI) unit member. All filenames are printed 
in a slanted font, rotonef.exe for example. Analytical symbols an' printed with a math italic 
font , .r for example. 
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Chapter 2 

ACOUSTICS REVIEW 

I his chapter gives a brio! discussion of acoustics concepts, to provide the non-acoustician 
with a quick reference for ROTONET work. Additional references for general acoustics 
include Pierce' [(>] and Kinsler et al [7], Various noise' metrics used in rotorcraft acoustics 
are' discussed by Pearsons Benmett [8]. Signal processing considerations arc' presented 
by Hardin [9] and Bendat and Piersol [10]. An overvienv of rot.or noise' is pre'sented bv 
Hubbard [11] and Brentner and Parassat [12]. The theoretical me'thods actually employe'd 
in IK) I ONE I are discussed in rele'rences [1], [2]. [d], [1]. and [r>]. Discussions of rotorcraft - 
specific noise' phenomena and issue's are' given in Chapter T 

2.1 Concepts 

hoi practical purpose's, unsteady pressure disturbance's in a compressible fluid is a clef- 
i i i i t ion of sound, hor aeroacoustics. the fluid is air. Noise' may be' dedined as annoying or 
undesirable sound. Sound wave's traveling through a fluid cause' fluctuations in pressure, 
density, and particle 1 velocity. I lie fiedd of acroacoust ics is that subset of acoustics which has 
aerodynamic phenomena as the noise' source's. Some' example's of aeroacoust ic noise' source's 
are': jed c'nginc's, propedlers. lielieopter rotors, and sonic booms. 1 lie' noise' of conce'rn may 
be' ewtorior t.o the 1 aircraft, such as community response* to airport noise', or inle'rnal to the 
aircra.fi. such as propelh'r noise 1 experiene’C'd by a passenge'r in a t urboprop commute'!' airline 1 !*. 
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Acoustic problems arc typically described in terms oi the “source . the “path , and the “ob- 
server." The source* is defined as the aerodynamic body or fluid parcel which originates the 
outgoing pressure waves. Because* the speed of sound is finite, a pressure wave receivc'd at 
an eebserver location at time / must have originated at the source at an earlier time r. 1 his 
earlier time r is commonly denoted the* “emission’ or “retarded time, while the time / is 
commonly called the "reception" lime. The path is the trajectory of the part of a pressure 
wave which reaches the observer. A differential surface element oi the pressure wave arriving 
at an observer has an incident direction, which in a two-dimensional problem is described 
by a reception angle. Likewise, that differential surla.ee element ol tin* pressure wave has an 
emission angle at the source. Observer may mean an actual person, a pressure transducer 
which is part of a data acquisition system, or a hypothetical transducer in an analytical 
representation of the pressure field. In any case, the observer is considered to exist at a 
single* point. 

Idle behavior of sound waves in a fluid can be modeled using the mass, momentum, and 
energy conservation equations lor a compressible continuous fluid, for low intensity sound, 
which is usually the case far from the noise source, the disturbances can be* considered 
mathematically small and allow assumptions leading to a linear governing equation. 1 he 
acoustic farfield is defined to be the region where the acoustic pressure gradient and the 
acoust ic velocity are parallel and in phase*. One form of the linear wave* equation is 


V-p = 


1 c) 2 p 

~ 2 W 2 


where* p is the aeemstic pressure*. / is time, c is the sound spe*e*el, anel X 72 is the Laplacian 
e>pe*rator. The time* variation of the* pressure field of an aircraft as me*asure*d on the* ground 
is elue* to motion of constant -strength sources relative to the* fixe*el ground ejbserver. as well 
as time variation of the seanves on the* aircraft . Ten* reeteere ralt . the* relative se)iirce/observe*r 
motion has t we> component s: the relative* velocity betwe*en the* airnali anel a ground observe*r. 
anel the relative veleieity betwee*n a pressure* source on a rot e>r blaele* anel an obse*rve*r in a 
non rotating frame*. 
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Various lorins of t. ho linear governing equation for sound have' been studied analytically 
for quite a long time. Since* the behavior of the linear equation produces ideal mathematical 
wave* phenomena, most acoustics work is couched in the mathematics and terminology of 
wave* mechanics. There are three common types of wave's considered in linear acoust ics. Plane* 
waves have flat wave* fronts; a classic plane wave fielel occurs in front of a flat piston mewing in 
a cylindrical tube. ( ylindrical waves are caused by phenomena like* an infinite-lengt h ewlinder 
expanding and contracting radially. Two-dimensional analytie-aJ problems such as oscillating 
airfoil flows involve cylindrical waves. Spherical waves are t hroe'-dimensional and are caused 
by source's like expanding and contracting sphere's, infinitesimal point pressure secures and 
the like*. Phenomena involves! in rotorcraft system noise studies are predominantly three- 
dimensional in nature. 

Acousties problems may also be classified into fre*e-field and reverbe*rant cast's. A free- 
hold ease is one in which waves emanating from a source* propagate* away from the source 
into the surrounding space without striking any solid object. Reverberant fields are those* in 
which senind wave's reflect off other solid bodies in the field. Examples of reverberant field 
problems are* aircraft interiors, noise reflection off the ground, an el refl(*e*t ion from buildings. 
Solid surfaces which reflect sound wave's may have either finite* impedance (they absorb 
senne* acoustic energy from the incident wave) e>r infinite impedance (all acoustic energy is 
reflevted). In a reverberant problem, the sound wave's can refle*ct e>ff of a solid surface*, or 
eliffract around an obstacle. 

Waves traveling in a fluiel e*an also be subject to refraction (a “bending" e>f the* wave's due 
to local change's in the fluiel properties such as sounel spe*e*d). scatte'ring (a disorganization 
of the* wave's due* to irre'gular fluiel motions or suspendeel sol iel part ides), and alternation 
el lie* to viscosity and e>t -her factors. In aircraft noise probleuns. refraction is fre'epie'ntly pre*se*nt 
elue t e> vertical gradients in air |>rope*rties (i.e. temperature* or wind ve'loe-ity) across the* 
altitude's be*t we*e*n the* aircraft and the ground. Since the* Kart Ids atmosph(*re* usually has 
some* t urbule'tice* pre*se*nt. scattering can be* a p]>r<*e ia ble* for aircraft flyove'r noise*. Attenuation 
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is always present, since air is a viscous fluid. All three of these phenomena (refraction, 
turbulent scattering, and attenuation) tend to become important for wave propagation over 
long distances, particularly high altitude flyovers. 

One of the difficulties ol aeroacoustics is that the sound source, usually an aircraft, 
is moving relative to a fixed observer. Most acoustic theory and technique is oriented to 
“stationary*' problems, for which the descriptive statistics of the acoust ic pressure (i.e. mean 
and variance) do not vary with time. For situations with relative motion between the source 
and observer, the problem is usually nonstationarv. An additional difficulty in aeroacoustics 
is that real-world problems have observers fixed to the ground, an aircraft flying on some 
(light profile, and some amount of mean and turbulent ambient wind. Thus, there are 
three 1 coordinate 1 systems involved rather than the 1 usual single 1 system, which complicates 
the formulation of modeling equations, and limits our ability to use physical intuition to 
understand the behavior of sound waves. Fortunately, many flyover situations consist of 
an aircraft flying at constant speed and altitude, in little or no wind, allowing a simpler 
formulation of the problem. 

The usual intention of acoustics work is slanted towards the end effect on humans. 
Undesirable sound is studied to determine the physiological and psychological effect caused 
bv a given source. There are non-lniman effects ol sound, as well. Acoustic fatigue of 
structures is one such effect, as well as acoustic modification of aerodynamic boundary layers. 
Though aircraft vary tremendously in size (thus their acoustic signatures vary substantially 
in magnitude and frequency), humans have relatively fixed length and frequency scales, so 
the end result of acoustic work is usually phrased in a dimensional or quasi-dimensional 
form. 

2.2 Units 

Sound, being a pressure disturbance 1 , has the units of force per unit area. Some units 
and conversions common in aeouslics are: 
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1 //bar 
I psi 

I (Ivno/an 

I Pa 


= 0.1 N/ni 2 
= 0894 N/m 
= 10-' N /in 
- I N /in 2 


Because t li<* dimensional magnitude of sound waves varies dramatically, acoustic pressure is 
often logarithmically scaled using the decibel (dB). Sound Pressure Level (SPL) is defined 
as 

SIM. = ailog.J^j dli (2.2) 

= “0 1og l0 p + 91 (2.4) 

when' is sot to 1 he commonly used reference pressure ‘20 fiPa. Often, the acoustician is 
only interested in magnitudes, so the mean square pressure, { p 2 ). is used: then 

SPL= 101og lo ^j dB (2.4) 


2.3 Frequency-Domain Representation 

Duo to the wave nature of sound and the many analytical and experimental tools appli- 
cable to signals expressed in sine series representations, pressure time histories an 1 typically 
Courier transformed to obtain amplitude and phase as functions of frequency. Plots of these 
functions are called spectra, with the amplitude vs, frequency spectrum used frequently in 
aeroacoust ics. Since the acoustic pressure time series are analyzed on digital computers, 
t -lif' analog signal from a microphone is digitized, converted from units of voltage to units of 
pressure, then stored for subsequent Fourier analysis. Note that in the present work, “band** 
is used to describe the width of the span of frequencies corresponding to a single Fourier 
coefficient . The term "bandwidth’* is often used to describe the width of the total spec- 
trum. which can lead to confusion. The frequency resolution of a spectrum is measured bv 
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the width of the frequency bands which is dependent on the digitization rate. Narrowband 
spectra are spectra which have bands that are constant in width (i.e. range of frequencies) 
and narrower than 1 /•? octave bands. Frequency resolution lor most rotorcralt narrowband 
spectral analysis of full-scale acoust ic data ranges from I to 10 Hz. Broadband analysis typi- 
cally uses —octave bands, for which the center frequencies (and thus the width ol individual 
bands) are logarithmically spaced. The advent of high-speed analog/digital converters and 
hast Fourier Transforms ( FFTs ) has made narrowband analysis a common tool. 

One analytic form for the Fourier transforms are 



though the reader will frequently find different ways ol handling the 2 tt factor in t he literature. 
As long as the Fourier integral identity holds, 

the form is correct, but care must be taken when comparing Fourier coefficients ironi two 
separate analyses. Since the p(i) time series is usually in digital form, tor both predictions 
and measurements, the Fourier analysis must be done in a discrete* lorm. I he general Discrete 
Fourier Transform (1)1 1 ) has the form 

tv-) = ^ £ cm 

ti — — oo 

However, since we are limited to analysing a finite lengt h time series, the 1)11 is restated as 

A/ N 

TV) = — £/>(» A/ )c-'^' (2-9) 

1 7T , 

n=l 

which can he applied to a finite time series of arbitrary length VAl. Since this transform 
is computationally intensive when applied to many time series with large \. various hast 
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Fourier J raiislorni ( I' I* T) algorit hms have been developed ior use in computer programs. 
1 liese ITT methods are very efficient, but most require that A be a power of 2. There are 
many issues requiring treatment when performing Fourier analysis, such as aliasing, leakage, 
and windowing; the reader is referred to Hardin [9] or Bendat and Piersol [10], 

When rotorcraft noise is measured and plotted in an amplitude spectrum, certain indi- 
vidual frequencies have magnitudes well above the rest of the spectrum. These frequencies, 
called tones, are generally considered to originate from deterministic, discrete phenomena. 
Other regions in the* frequency spectrum have relat ively equal magnitude across a slioii range 
oi frequencies. The content, of the acoustic signal creating these Hat spectral regions is called 
broadband noise. Broadband noise tends to come from non-deterministic sources, usually 
turbulence in sonic form: the viscous boundary layer on the blade, aerodynamic response to 
atmospheric turbulence, turbulence in the rotor wake, and turbulence in the engine plume 
are' all cause's of broadband noise'. Broadband noise' can be' a.de'ejuate'ly describes! in the 
frequency domain with ^-octave bands, while' tone' noise' ie(|uiies narrowband analysis to 
properly capture' the 1 height and shape e)f the' tone' in the' amplitude* spe'ctrum. 

Human hearing is generally consielered to span fre'que'ncie's from 2 Ilz to 20 kHz. He>w- 
e'vc'r. for most aeroacoustic applications, the upper limit is induced to 10 kHz. Since 1 much 
signal processing is performed with ^-octave band analysis, the uppe'r limit is e'xtende'el to 
1 1.2 kHz. which is the 1 upper limit of the 10 kHz ^-octave band. For jet noise, the 1 lower limit 
is ge*ue*rallv set to 50 Hz. but rotorcraft often have dominant tonal eontent bedenv 50 Hz which 
should not bo overlooked. Note that, since 1 1 lie* end result is human perception of the 1 noise, 
scaling effects heroine' important. If one is interested in estimating the 1 noise of a full-scale 1 
rotorcraft by measuring the* noise 1 from a sub-scale' model, t lie free|iie'ncie's of inteirsl must 
be sliiite'd. For example, if a one'-fourth scale' model operates at constant tip Mach number, 
t lie'll the* rotor rpm would incre'ase by a factor of four, and likewise' the* characteristic noise* 
frequencies. T hem if the noise* of the full-scale aircraft- is re<|uired up to 10 kHz. then the' ~ 
scale 1 model noise* must be* measured te) 10 kHz. 
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To look at a frequency-integrated sound level, the Overall Sound Pressure Level (OASPL) 
is defined by 

OASPL = l()log 10 £ dB (2.10) 

" \P ref J 

where (p 2 ) n is t he mean square pressure of the //-ill frequency hand. In an at tempt to quantify 
the frequency response of the human ear, particularly increased sensitivity to high ( 1-2 kHz) 
frequencies, various frequency weighting schemes have been developed. Two common ones 
in use for aeroacoust ics are A-weighting and PNL weighting. I liese schemes multiply the 
magnitude of each frequency hand by some weighting factor, then sum the products as 
for OASPL. An additional corn'd ion to the PNL scheme is the PNLT (Perceived Noise 
Level, Tone-corrected). This metric weights distinct tones and adds their contribution to 
the PNL- weighted broadband content of the spectra. This additional step is performed to 
quantitative! v model the additional annoyance that is associated with tonal sources, for 
aircraft flvover noise, tin* Effective 1 Perceived Noise' Level (KPNL) is used as a certification 
met ric [13] as well for est imation of community response to noise levels. KPNL is obi aim'd 
by a time integration of PNLT during a flyover. 


2.4 Typical Noise Levels 

The' following table presents some' typical noise levels in decibels and pounds per square 
inch to illustrate the effect of log scaling and give' a reference for aircraft noise relative to 
ot her noise' sources; 


noise' 

dB A 

hearing threshold 

0 

whisper 

20 

quiet neighborhood 

10 

normal speech 

00 

busv office 

so 
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heavy traffic 100 

discot heque 120 

where (IB A means the sound pressure* levels are A-weighled. Note that a helicopter flying 
overhead at 500 feet altitude may have a noise level of 80-100 dBA. Even at noise levels of 
110 dBA. which is a rough estimate* ol the pain threshold, the root -mean-square acoustic 
pressure is only 200 Pa. which is a very small variation about the 101. 225 Pa static pressure. 
Noise levels of 100 dB. which is loud enough to make conversation difficult, have a root -mean- 
square acoustic pressure of only 2 Pa. This difference in magnitude between the mean and 
acoustic pressure lias severe ramifications when trying to directly compute acoustics from 
governing fluid dynamics equations, since numerical errors can be larger than the acoustic 
content of the calculations. 

2.5 Air Properties 

Some properties of air at sea level, standard conditions are listed below. Properties of 
the V. S. standard atmosphere for different conditions have been tabulated in English and 
metric units for a series of altitudes [1 I], 

pressure 10 1225 Pa 

density 1.225 kg/m* 

sound speed 310.29 m/s 

I sing A = j to relate wavelength A to sound speed c and frequency /'. the following list, 
illustrates wavelengths in the frequencies of inlerest for a sound speed for air of 310 meters 
per second: 
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f (Hz) 

A (meters) 

20 

17 

200 

1.7 

2000 

0.17 

20000 

0.017 


since helicopter tone noise occurs roughly between 15 and 1)00 Hz. the lowest frequencies 
have a wavelength larger than most helicopters, while the highest are* smaller than the rotor 
blade thickness. 
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Chapter 3 


HELICOPTER NOISE 


Helicopter aeroacoustics has become a distinct field in recent years, due to the increased 
interest in FA A certification and military survivability issues. The aims of the workers in 
this field are to understand the physical mechanisms which cause' noise, to measure and/or 
predict, the noise level of a particular aircraft for comparison purposes, and to reduce tin' 
noise level of future designs. The tremendous expense and difficulty in getting detailed, 
accurate experimental data has hampered efforts on all three fronts. Recent experimental 
programs have sparked substantial progress in understanding the aeroacoust ic processes, and 
have* illustrated the* need for improved analysis and design tools. This chapter discusses tin* 
various noise source's on rotorcraft and the current methods lor predict ing them. 

3.1 Overview 

The noise* signat ure* generate'd by rotorcraft is quite complex, and has manv difficult - 
to-moelol source's. Aerodynamic noise' is gene'ratod from main rotors, tail rotors, fuse'lage 
boundary layers, and engine's. While* aeroacoust ics may be thought of conceptually as line* 
resolution of unsteady pressure's in an aerodynamic flowfield, it is pragmatically difficult to 
modc'l the* endin' flow re'gion, which includes the* helicopter, ground plane', and observers, 
f bus, pre'dietion efforts to date* have* klent ifie'd the* more important noise* me*chanisms and 
applied focuse'd approache's on more* tradable* models. 
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Some aerodynamic phenomena on rotorcraft are: unsteady self-interacting vortical ro- 
tor wakes, transonic flow around the blade tips in high-speed forward flight, rotor/fuselage 
interaction, rotor/rotor interaction, blade boundary layer separation, continually varying 
blade pitch due to pilot -commanded control inputs, rotor wake/engine plume interaction, 
and atmospheric turbulence. Likewise, computing propagation of the relatively weak acous- 
tic waves using nonlinear computational methods has been very difficult. The solution of 
some' form of the governing equations for unsteady compressible fluid dynamics lor the re- 
gion enclosing the vehicle and observers is a task which will not become practical for some 
time, due to limitations of current algorithms for solving the governing equations as well as 
limitations of available computer hardware. 

The Lighthill acoustic analogy is used to separate the aerodynamics from the acous- 
tics [15, 1(>], In this way. the aerodynamics can be solved using existing methods without 
requiring the resolution needed to capture acoustic waves. I he aerodynamic modeling typi- 
cally used reduces the problem to that of an isolated rotor, since the interaction with other 
bodies and rotors is an extremely difficult problem, and t he isolated rolor problem is the pri- 
mary contributor to the total noise signature. Using the aerodynamic loading on the rotor, 
the generation and propagation of sound is then predicted with several acoustic theories [12]. 
The noise has contributions from both rotor tone noise and rotor broadband noise models. 
A recent state-of-the-art review article on helicopter rotor noise is given by Brentner and 
Farassat [12]. Another overview of rotor noise is presented by Hubbard [11]. 

3.2 Tone Noise Sources 

A substantial understanding of rotor tone noise has been developed using the I'towcs 
VVilliams-llawkings (FW'-Il) equation [17]. One form of the FW-II is 

hw - v ‘" = iJi 1 1 - it, it I ^ 1 - 




[T„ n ,n(f )} (3.i) 
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I here are several forms of the hYY-II equation to he found in the literature [18, 19, 20. 21], 
1 his equation, which is a rearrangement of the fluid dynamic conservation equations, can be 
reduced to a linear three-dimensional wave equation. 'The lei t hand side is a wave operator, 
and the right hand side has forcing functions due to the kinematics and aerodynamics of 
the problem. The kinematic forcing term can be construed to be the effect of a finite- 
thickness blade displacing fluid as it passes through the air, thus it is called the “thickness 
term. Likewise', the aerodynamic term is called the “loading term" since' it is due to the 
aerodynamic force (or load) on the fluid by the black'. This nomenclature 1 can be misleading, 
insofar as a nonlilting ease still has loading noise, due to the surface 1 pressure 1 and slic'ar 
distributions on a blade 1 section at its angle of zero lift. 

A Greens function solution to the FYY-II eepiation has been developed bv F. haras- 
sat [18]. I lie 1 solution is known as “formulation 1A and is given bene, 

p{.vJ) = p L (.rJ) + pr (.vJ) ( 3 . 2 ) 


Jr=o 


where the 1 thickness noise 1 pr(xj) is 
l7r/> r (.r. /) = 

and the 1 loading noise 1 is 

I xpi.{r.l) = ~ I 

C JS 


Pnt 


'ii "I" cM r — c.\l 2 j 


r 2 ( \ - Mr) 


(IS 




hr, 


S=1> !_/■(! - Mr) 


<IS + / 
Js = 0 


rd - Mrf. 


(IS 


+ 


c JS-o 


l r (I’MjVj + <M r — C.\l -'j 


Mr) 


(IS 


(••M) 


rrt 


I he thickness noise, i.e. noise predicted by the thickness term, of a rotor is predomi- 
nantly in-plane with the rotor and is highly directional. It is affected by the rotor tip speed, 
advance ratio, number o( blades, and airfoil shape. Since it is based only on kinematics and 
geometry of the blade, prediction methods have been successful for this contributor to tom* 


noise. 
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Loading noise effects are strongest out of plane with the rotor, i.e. above and below 
the rotor disk. Loading noise is affected by thrust coefficient, advance ratio, number of 
blades, blade twist , blade camber, tip speed and rotor t rim. The aerodynamic forcing to the 
FVV-H equation is typically obtained for an isolated rotor encountering level steady flight 
conditions. Even for this ideal situation, only recently have the rotor aerodynamicist s tools 
been adequate for use in predicting torn' noise. Some preliminary work has been done on 
rotor/airframe and rotor/rotor interactions, but these techniques are still exploratory in 
nature. 

The effects of compressibility on the acoustic field of a rotor are twofold: the noise 
due to shocks, and the alteration of wave propagation due to varying air density near the 
blade. They may be lumped into the High-Speed Impulsive noise (HS1) category. Both of 
these phenomena are modeled by the quadrupole terms of the FW-H equation and are not 
commonly modeled. Research is underway to determine how to incorporate their effects into 
FW-H equation solvers. Additionally, a technique called the Kirchhoff method is being used 
to predict far field rotor noise by direct- computation of the pressure field around a rotor blade 
using Computational Fluid Dynamics (CFD) codes for a small region surrounding the blade, 
and a linear wave equation method for wave travel beyond the the (TI) solution region. In 
this wav. nonlinear methods (i.e. the inherent nonlinearity of the (TI) governing equations) 
are used to directly compute the pressure wave behavior in regions with nonlinearity, and 
an eflicienl linear wave equation where the problem is predominantly linear [22. 23]. 

Blade- Vortex Interaction (BYI) is an impulsive aerodynamic phenomenon where a well- 
defined tip vortex passes near a rotor blade, and induces sharp spatial and temporal gradi- 
ents in the blade surface pressure' distribution. These surface pressure fluctuations generate 
high I v directional acoustic impulses which can greatly increase the peak-to-peak pressure 
experienced in the farfield. A main rotor can have B\ Is with its own wake in high speed 
level steady flight, and also in forward flight descents at intermediate speeds. In high speed 
forward flight, the noise problem can be magnified by the presence of shock waves on the 
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blade t.i|> during the advancing quarter o( each revolution. The vortex interaction with the 
transonic How around the blade tip can generate very strong impulses. Tail rotors encoun- 
tering main rotor vortices also have BVIs. The propagated noise for both causes is highly 
directional and extremely unpleasant. 

f he tones generated by a helicopter rotor occur at frequencies which are integer mult iples 
ol the "Blade Passage Frequency’* or BPF. which equals the number of blades multiplied 
by the rotor angular speed in Hz. The BPF for a McDonnell Douglas MD-500 helicopter is 
aboul 11 Hz, which is the vehicle used lor the example cases discussed later in this document. 
For flight conditions in which there is no BYT or shock noise, the tonal part ol the spectrum 
is appreciable lor the first >-10 harmonics ol the BPF. Impulsive' phenomena such as BY I 
and IIS1 increase the tonal content beyond 30 harmonics [21]. 

3.3 Broadband Noise Sources 

There are many mechanisms which contribute broadband noise to rotorcraft acoustic 
signatures. Turbulence ingestion noise is generated by the fluctuation on rotor blade loading 
caused by ingestion of atmospheric turbulence and the resulting azimuthal change in blade 
effective angles-ol-attack. Kngine broadband noise arises from the turbulent mixing of the 
exhaust stream with the freestream air, the combustion noise inside* t lie* engine, and com- 
pressor noise. Another broadband noise source is Blade- Wake Interaction (BVYI) [25]. This 
is caused by turbulence in the blade's shed wake interacting with the black' aerodynamics. 
Very little research has been done on this topic. Turbulent boundary layer trailing edge 
noise 1 is caused by the passage of turbulent, eddies in the boundary layer passing bv the 
blade trailing edge. Laminar boundary layer trailing edge noise is due to vortex shedding 
at the trailing edge' from the laminar boundary layer. Separation-stall noise can lx* caused 
by both mild trailing edge separation regions as well as massively separated regions. Very 
little' is known about this noise source'. Vortex formation noise is due to the rapid creation 
ol turbulent flow in the formation region of a rotor tip vortex. Trailing e*dgr blunlness noise 
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is caused by vortex shedding horn a blunt trailing edge. These noise sources are discussed 
in more detail by Brooks it al [25. 26], 

While current technology rotors have amplitude spectra dominated by low frequency 
tones. freq uei l cy- integrated noise metrics show that broadband noise is as important as tone 
noise. In fact, broadband noise tends to dominate the total levels computed using frequency- 
weighting (dBA, for example). However, prediction of broadband noise is a difficult task, 
since turbulence and/or chaotic phenomenon; thus, current technology for the prediction oi 
broadband noise is limited to combinations of semi-empirical and analytical methods. 

3.4 Propagation 

As sound waves pass through the atmosphere, then" are losses in signal strength due to 
several causes. The dominant effect is "spherical spreading", which is due to the increase* 
in area of a spherical wave as the wave front travels away from the source". Note that 
there is no loss of total acoustic energy, just a reduction in energy per unit area per unit 
time. Another effect is the loss ol acoustic energy by heating the fluid and other internal 
energy loss mechanisms. These losses arc* generally denoted "atmospheric absorption. Tin* 
effect of a wave striking the ground can be large. The wave is partly absorbed and partly 
reflected, with a reflected signal arriving at an observer some lime after the original wave. 
Atmospheric turbulence can also scatter sound waves. Lastly, temperature, humidity, and 
wind velocity gradients can refract, the sound waves. All ol these ("fleets attenuate the original 
noise generated by the aircraft. 

In the human-audible range, the attenuation is frequency-dependent, with greater at- 
tenuation at higher frequencies for most problems. Note that the presence of the ground 
vastly alters the propagation due to diffraction and thermal gradients. Propagation in free 
air is a much simpler problem. (Jenerally. 100 meters is a typical distance beyond which 
atmospheric effects are considerable*. 
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Some means of combining propagation mode'ls and source' models is needed for rotorcraft 
flyover noise prediction. 1 he method used in ROIONLJ is to use' the 1 noise' source' models 
to predict the* noise' in the' moving aircraft frame', with the' only propagation effect In'ing 
the' spherical spreading behavior implicit in the governing expiations, and use' separate' proj)- 
agation models to predicl the' attenuation, scattering, and ground reflection. The' transfer 
of information between the source* models and atmosphe'ric propagation models occurs at a 
hypothetical "‘sound bubble' , a spherical surlace some' 5-10 rot, or radii away from the* rotor 
hub. 

3.5 Noise Metrics 


Major driveTs lor developing rotorcraft noise' pre'diction capability are' the eommunitv 
noise' and cert ifieation problems. Two metrics are' commonly use'd to quantify the relative 
human response* to aircraft flyovers. The* EPNL metrie' is a PNL-based (piantitv for a single* 
flight event, while the Day-Night Level. DNL (sonietime-'s denoted by Lu\ ) is an A-we'ighting 
base'd metric for an assumed seric\s of flights during both daylight and nighttime hours of a 
21-hour day [8]. NYithe'rof the'se' sen'm to properly asse'ss the' annoyance' of a low-fre'ejuencv 
dominated vehicle such as the' UII-1. Research is ongoing in the' psychoacoustics fie'ld to 
dewedop better models for the* relative human response to tone' noise' and particularly, highly 
impulsive* noise'. 
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Chapter 4 


ROTONET EXECUTION 


ROIONKT is a large Fortran program which reads in the flight conditions and a user- 
created helicopter database, and execute's a user-specified sequence ol subrout ines to predict 
the noise' received at ground observer locations. This computer program, developed bv the 
Acoustics Division of the NASA Langley Research (’enter during the 1980's, (‘an be installed 
on several computer types and operating systems. The complexity of the problem and the 
diversity ol models involved in utilizing ROIONKI for prediction and design work requites 
a substantial commitment of time to learning the ANOPP executive and the mechanics of 
the various ROTONKT modules. This and oilier publications [1. 2, .*{, 1] provide the user 
with information needed to run ROTONKT. This chapter discusses the way execution paths 
are configured in order to make various kinds of noise prediction with ROTONKT. 

4.1 Basic Concepts 

R010XK1 actually consists of a subset of the ANOPP Aircraft NOise Prediction Pro- 
gram. ( ertain user-invoked modules are helicopter-specific, and, when given the proper 
inputs and execution sequence, product' the desired outputs. Depending on ihe computer it 
resides on. ANOPP may be run either interactive- or batch-mode. For either, a procedure 
or script file is required to assign the input and output files to logical I/O units, as well as 
run the execut able file. 
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The ROTONET program is modular. To execute certain functions. Functional Moduli s 
are invoked in the input deck. Each ol these modules is a sel of subroutines which do the 
various calculations. Required data are fetched, operated on by a specific algorithm, and 
results stored in ROTONET s run-time database as well as printed out to the standard 
output file. 

For anv given ROTONET run. an input file (also calhnl an input deck) is needed to 
provide the code with input data and an execution sequence for the chosen modules. Example 
decks are illustrated in later sections ol this chapter. 

The many modules in the ROTONET group can be grouped in the following manner: 

• atmospheric conditions source/observer geometry 

• blade* properties 

• rotor performance* 

• rotor tone noise* 

• tone noise propagat ion 

• rotor broadband noise 

• engine' broadband noise 

• broadband noise propagation 

• noise met ric calculation 


4.2 Functional Module Groups 

Table 1-1 lists the modules used in typical R.OTONET analysis, as well as a brief de- 
scription of their functions: the* modules are* grouped by similar and/or related functions, 
and are listed in order of execut ion for most analyses. Figure 1.1 presents a flowchart of the 
entire ROTONET system. 
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Module Purpose 

ABS atmospheric absorption coefficients 

S VO steady level flight trajectory 

(iKO source/observer geometry 

APM atmospheric property parameters 

ATM atmospheric properties 

IPS reads and interpolates blade shape 

1BA inviscid aero forces and moments 

IBL viscous forces 

LHP rotor trim and aerodynamics 

HYY(» rigid wake geometry 

H IN wake-induced velocity 

H HI) rigid-blade flapping response 

HU) rotor airloads 

LHX compact-chord wise rotor tone noise 

HBN rotor broadband noise 

HTN full-surface rotor tone noise 

ABL atmospheric turbulence properties 

S I L streamline distortion tensor 

LIN turbulence ingestion noise 

MSN combines source noise for two rotors 

PHO propagates rotor broadband noise 

PHI propagates rotor tone noise 

LKY (requeue v-integrated noise metrics 

EFF EPNL for <*ac 1 j ground observer 


I able 1.1: ROTO NET Modulus listed bv fundic 



UOTONET Primer 



Figure 1.1: KOTONEI System Flowchart 






























HOTONET Pi inter 


4-5 

The acoustic properties anti source/observer geometry generate inputs to the propaga- 
tion and metric calculation modules. The sequence of executed modules is ATM-ABS-SFO- 

cfo. 

The blade properties modules generate inputs for the rotor performance and rotor broad- 
band noise modules. There are two sets of modules which can be used; either RBS-RBA 
BLM . or I BS-IB A - 1 BR. I he former are the original modules, and t he lat ter art* more recently 
developed modules. The improved modules should be used whenever possible, and the orig- 
inal modules will not be discussed further in this document. 

The rotor performance modules create airloads tables for t he rotor tone noise and broad- 
band noise modules. For simple blade element /momentum theory calculations, the LRP 
module is executed, for higher harmonic loads, the RLI) module is executed, in conjunction 
with HR I). RIN. and RYYO in an iteration loop. 1 lie rotor tone noise prediction generates 
pressure time histories or spectra lor direct use or lor input to propagation modules. If the 
LRP module has been used, the LRN module must be exercised; if rotor loads were generated 
with the RLI) system, the Rl N moduli' must be used. The rotor broadband moduli' RUN 
likewise generates spectra. 4 hose may be used directly or as input to the broadband propa- 
gation module. Rotor tone noise propagation is computed with PR L. Its calculations may be 
used directly (narrowband spectra) or input to the noise metrics modules. Broadband noise 
propagation is calculated with the PRO module. As with PILL. PRO's calculations may 
be used directly (^-octave band spectra) or input to the noise metrics modules. The noise 
metric modules are called with the sequence LEY-KFF. Various integrated noise metrics at 
a ground observer are produced from the propagated noise generated by the modules PR I 
and PRO. 

I he various Functional Modules are discussed in detail in Chapter 5. Their inputs and 
outputs are discussed in Chapter 6. 
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4.3 Input Deck Structure 

When ( real ins a KOTONKT input, deck, the following choices must be made; 

• what observer locations and reference frame are desired 

• what noise outputs are needed (EPNL, --octave' band spectra. 

narrowband spectra, time-histories) 

• what noise sources are to be modeled (main rotor tom' noise, main 

rotor broadband noise, engine broadband noise, etc.) 

• how rotor aerodynamic loading is to be predicted (simple momentum. 

generalized wake + rigid rotor dynamics, experiment) 

By working backward, the required modules arc' found, which ihen determines the order of 
execution via the needed output-to-input relationships of the module's. Also, the external 
user inputs are defined at this point. 

There arc' several input deck template's which illustrate workable ROTONET models. 
Often, the fastest way to create' a new input clock is to reconfigure one of the templates for 
the problem at hand. However, the flexibility of the executive structure is such that a great 
many more possible' configurations may be used to accomplish a new task. 

4.4 Inputs and Outputs 

There' arc' two basic types of inputs and outputs; parameters and user database members. 
Both may be used as input, out put , or both. Parameters are single valued or one-dimensional 
arravs. while' members arc 1 arrays, frequently two- and thre'o-dime'nsional. Note* that param- 
eters and members are global during the' run; the' information is in a common memory area 
and all functional module's and data manipulation commands have access to the* value' ol a 
database 1 element . 

Paramc't ers are' integer, real, logical, or character data arrays oi one* or more* elements. 
Most parameters arc' single* element real numbers. Some* are* onc'-dimensional arrays. Most 
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parameters have default values as well as tipper and lower limits. 

Database members may be unit, members, table members or procedure members. They 
can be either “run-life" or part of a database' stored on disk. Note that. AN01T has a 
permanent database containing all data ANOPP provides for itself. A database can consist 
ol several members, each of which is composed of records. Kach record can consist of one or 
more' elements, hlements can be character strings, integer, real, or logical data. 

P nit members are database members ior data which are' used as-is, and not interpolated 
horn. An example' is the GRID(XIl) member which represents a spanwise blade* grid distri- 
bution. 1 able' members are database members which are' designed to be interpolated from. 
Wlie'ii a table' is created, the use*r specifiers the type' of interpolation useel. An ('xaniph* is the* 
EMF(PRES) table', which consists oi a blade* surlace pressure* clist ribiit ion on a give*n grid. Ihe* 
pressure's are inter polated to obt ain the* pressure' at t he* elesired comput at ional blade* grid. 
Proce'dnre members are* portions of ANOPP input de'cks which may be executed during part 
ol an ANOPP run. One* may think ol them as a "macro , or se*t of commands which are' 
free|uently used in the* same' sequence*. 

A thorough discussion of the ANOPP input syntax is given in reference [2]. Definitions 
of input and output variable's are given in "The NASA Aircraft Noise* Pre'diction Program 
Improverl Pro])e , lIe*r Analysis System [-5], "Aircraft Noise* Pre'diction Program The*oret ical 
Manual volume* 1 [l]. An Introduction to High Spe*e*d Aircraft. Noise* Predict ion" [5]. and 
"Aircraft Noise* Pre'diction Program User's Manual" [2]. Another excellent refere*nce* for input 
anel output variable's can be* formed by making a printout in alphabetical order of the* “help 
file's ( V ax ve*rsions) or "man page's (Unix versions) lor all ol the* module's use*d in typical 
KOTONKT predictions. 

4.5 ANOPP Control Statements 

ANOPP control statements are* available* to specify or modify input data, execute mod- 
ule's. manipulate' data arrays (eilhe*r in "unit me*mbeT or "table me'inber" form), and perform 
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general-purpose logic functions such as iteration loops and value tests. The following is 
list of the more commonly used control statements tor the ANOPP executive [2]; 


(General Executive 

ANOPP 

STARTCS 

TITLE 

ENDCS 

Parameter Controls 

PARAM 

EVALUATE 

UPLIST 


Library Controls 


LOAD 

UNLOAD 

DROP 

Dat abase Member Cont rols 
TABLE 
UPDATE 
PURGE 

Database List Controls 

CATALOG 

MEMLIST 

TABLIST 

Functional Module Controls 
EXECUTE 


starts ANOPP executive and sets some parameters 
initial executable control statement 
sets title string for standard output file 
(appears on printouts, can change during run) 
stops processing input deck and closes files 


used to assign a parameter to a fixed value 
computes the value to be assigned to a parameter 
causes printout of all exist ing parameters k values 

read in library from disk 

copy active library to a disk or tape file 

remove a library from the current job memory 

add a table member to t he database 

add or modify a unit member in the database' 

remove a specified unit from the database 

produce summary list ol current- database* units, 
or contents of a specific data unit 

give detailed information about a specific unit membe 
give* detailed info about a specific table* 

run t he specifier! mexlule 



ROTOSET Primer 


4-9 


I lie following section provides an example of actual deck structure and ANOPP Control 
Language syntax. The reader is directed to the Aircraft Noise Prediction Program User's 
Manual [2] for thorough details on control statement syntax. 


4.6 Sample Input Deck 


I he following input deck is used to illustrate the form of a "typical" ROTONKT deck. 
Note that it is a basic text file' (ASCII on most computers). Two library files, containing 
various unit and table members, are being read in. Thus, then' must have been one or 
more ROTONKT jobs previously run to create the library files. This feature is used to 
avoid running modules over and over for different flight conditions, but with the same blade 
geometry and properties information. Note that the library files arc binary, and cannot 
be read from a terminal. I he functional modules used here are described in the following 
chapter. 

ANOPP JECH0= . TRUE . JL0G=.TRUE. NLPPM=48 $ 

STARTCS $ 

SETSYS JECH0=. FALSE. $ 

$ 

$ example input deck for ROTONET Primer. 

$ 

$ note that on each input line, text after the dollar sign is ignored 
$ and is used for comments, note all caps for actual ANOPP commands. 

$ the dollar sign can be put in any column desired; Fortran style 
$ commenting is being used in this input deck. 

$ 

LOAD /BLADE/ $ a library file from a previous ROTONET run, 

$ with IBS/IBA/IBL results, is read into the 
$ database. BLADE. WRK is the file name on 
$ most machines, note all caps. 

$ 

LOAD /ATMSP/ $ a library file from a previous ROTONET run, 

$ with atmospheric and flight dynamics 
$ information. 

$ 

$ now define some parameters; let the rest 
$ default, note the spaces before and after 
$ the = signs. 
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$ 

PARAMETER R = 10.0 $ 

PARAMETER IMPROV = .TRUE. $ 

PARAMETER MACHRF =0.5 $ 

PARAMETER MZ = 0 . 1 $ 

PARAMETER NBLADE =4 $ 

PARAMETER WT = 0.006 $ 

$ 

EXECUTE LRP $ 

S 

EVALUATE BETA(l) = AOL $ 

EVALUATE BETA (2) = AIL $ 

EVALUATE BETA (3) = B1L $ 

$ 

$ 

EXECUTE LRN $ 

$ 

$ 

EXECUTE PRT $ 

$ 

$ 

UNLOAD /NOISE/ $ 

$ 

$ 

ENDCS $ 


blade radius, meters. 

use the IBS/IBA/IBL blade tables. 

hover tip Mach number. 

flight Mach number. 

number of blades . 

aircraft weight coefficient. 

compute the rotor airloads. 

convert LRP output parameters to LRN input 
parameters . 


compute the rotor tone noise in the a/c 
reference frame . 

propagate the tone noise to the ground 
observer . 

this writes a library file called NOISE. WRK 
for future use. 

this finishes the run, and closes the files. 


Any text following the ENDCS control statement is ignored. This feature allows stopping a 
job stream in the middle, by inserting an ENDCS statement in the middle of the input deck. 
This can be used to save CVV time, which is useful during debugging. 


4.7 Notes on ROTONET Use 


ROTONET is currently available on VAX/VMS and souk* l ’nix workstations. Included 
on the distribution tape are: an executable, help files (VAX/VMS) 01* man pages (l nix), cer- 
tain input decks called "templates” along with matching out put files (to verify the numbers 
generated on a user's particular machine), and an interactive' command file (VAX/VMS) 
or esh shell script (Unix) to run ROTONET. This command iile/shell script prompts the 
user for file names, creates a temporary working subdirectory, and runs the ROTONRI 
executable' with the specified input file's. 

The standard filename extensions on VAX/VMS systems are: (or the standard 

input deck (containing the 1 ANOPR control statements), "T/X\J for the' standard output 
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file (containing print or- format, ted stuff), Of 7" for the file containing array data in ASCII 
tabular form (from LRN when IPRINT=4, for example), and "*AYRK" for library files and 
temporary work files used by ANOPP for data management . In Unix, the filename extensions 
an 1 lower case except for library files, which must be completely upper case. 

Creation ol new input decks is a tedious process. The many “templates" and other 
demo decks provided with the ROTONET distribution tape provide' sufficient examples for 
handling most predictions. The process for creating new decks can be broken down into 
five' steps; eletermining the “end result” required, finding the sequence of Functional Moelule 
executes which will result in that end result, managing the data flow between modules, 
collecting and assembling the raw input data neeeleel to start the' various input elecks. anel 
debugging the elecks. The' first three steps are done with the help of this document. the' 
RO I ONE I theoretical manual [1], and the help code/man page's. Putting together the' raw 
elata needed tor moelules such as IBS can be' quite difficult. Preelict ie>ns are' epiite sensitive 
te) ejuantit. ies inclueling the geome'try of the' leading eelge' region of the' blade, the geometric 
twist rate* oi the' blade, the vehicle equivalent flat plate' elrag area, and so forth. 

There * are several steps to debugging ROTONET inpul decks. The first step is getting 
the decks to run through completely. IT the decks finish successfully, that is with no run-time 
errors, t he last two lines ol t he standard output file will he 

ENTERING ANOPP NORMAL TERMINATION PHASE 

ANOPP IS TERMINATING NORMALLY 

this only means that the ANOPP program finished without run-time errors. However. 
ANOPP may have had errors in reading the inputs or executing modules, in which case 
there will he error messages in the output file which can he found by string-searching in 
a texl editor lor ERROR. Since the standard output files can become very large, if is much 
quicker to use a text editor (in l.’}2 column mode) to find error messages. To cut down 
on the run time while debugging a particular stretch of inputs, a new ENDCS $ line can 
he inserted in the deck just after the problematic control statements. Thus, the rest of 
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the deck can remain intact and when the problem is resolved the temporary ENDCS $ line 
can be removed so that the rest of the deck is executed. ROTONET outputs can be very 
sensitive to several input quantities, including; the distribution of grid points on the blade 
(GRID(XIl) and GRID(XI2)), the number of time steps per rotor revolution (NTIME), the 
distribution of observers on the "sound bubble" (SFIELD(PHI) and SFIELD(THETA) ). and 
the number of ground observers ( OBSERV (COORD) ). It is highly recommended that t he user 
perform a study determining the necessary resolution by using increasingly fine grids until 
the noise predictions stabilize. There will be some variation with grid changes, but noise 
levels probably should not change more than 1 dB. If highly directional noise phenomena 
like thickness noise, high-speed impulsive noise, or Blade- Vorl ex Interaction become impor- 
tant, the observer arrays on the sound bubble and ground plane' should be reviewed. If the 
observer distributions are inadequate, the prediction accuracy will vary substantially with 
changing flight conditions. 

Due to ANOPP s origins in prediction for (TOL aircraft , the' LEV module expects the 
SFIELD(FREQ) array to contain the 21 standard —octave band (‘(Miter frequencies. The stan- 
dard center frequencies start at 50 Hz and end with 10,000 Hz. If the array does not match 
the standard 21 bands, an error message' is printed in the out put lile; ANOPP continues exe- 
cuting, however. Lor rotorcraft prediction, the noise 1 usually extends below the first standard 
--octave band centered at 50 Hz. Depending on the main rotor s blade* passage frequency 
(rotation rate in Hertz multiplied by the number of blades), the* SFIELD(FREQ) array may 
need to extend down several more* bands. LEV will continue on with the nonstandard array 
of band center frequencies, so users can ignore* the error message. For EPNL purpose's the 
21 standard bands (with center frequencies starting at 50 IIz and (Miding at 10 KHz) should 
be* use*d. Henveveu*, to correctly account for the low freepiency contend of main rotor noise, 
the' set of bands should be* e'xteuided down to the 12.5 Hz center Ireepiency. 

The* interpellation of the blade* shape* is a frequent source of problenns. If the distribution 
of airfoil ordinates is t oe> eourse* or too fine, e>r if the* GRID(XI2) array has excessively fine 
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spacing around t he loading edge (XI2 = 0.5), the IBS modulo lias difficulty interpolating the 
blade" shape. Also, it the trailing edge lias substantial camber, the transform to the 1 complex 
plane (discussed in section 10.2 oi Zorumski Weir [1]) may cause 1 the 1 angular coordinate 
to be'come non-monotonic, thus making interpolation impossible. The simplest fix to this 
problem is to remove 1 the camber from the trailing edge. 

One final comment pertains to deck maintonane*o and modification. Due to the large 1 
amount of input data, and the 1 difficulty in obtaining accurate inputs, it is recommended 
that te'xt be added to e'ach deck in ANOPP control statement “comment" form (see decks in 
Appendices ( \ 1), and K). Information to provide should include 1 ; which 1 type 1 . source 1 of input 
data. tein])ora.ry modifications for debugging or parametric studies, local point -of-cont act. 
file 1 name 1 , HO I ONI 4 , 1 system leaved, and purpose 1 of the 1 deek with list of modules called and 
libraries loaded. Archival versions of validated elecks should be 1 kept in a sale 1 plaee 1 se> that 
the v are 1 not overwritten in the 1 course 1 of a new project. 
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Chapter 5 


FUNCTIONAL MODULES 


Tlx- following pages contain a description of the functional modules in R0T0NET and 
tlie various coordinate systems involved. Commonly used ANOPP module's an' also briefly 
discussed. Detailed discussion of t he input/out put variables are presented in Chapters (> 
and 7. theory lor modules PRO. KFF. and LKV are discussed in Part I of the ANOPP 
Theoretical Manual [I], 'theory for RIPS, RRA, HLM. and PRT arc eliscusse'd in Part 3 
ol t fiat report . I heory lor the KO I ONK I modules is discussed in Pari 1. Usage' of the 
PHO. KFK, and LKV modules is discussed in the "Aircraft Noise Prediction Program User's 
Manual [2]. The 1 IBS, IBA. and IBL module s are discussed by Nguyen ['•]]. 


5.1 Coordinate Systems 

I here 1 arc' a variety of coordinate systems involved in a ROTONKT analvsis. Tracing 
from source to observer, there' are the; blade, shaft, tip-path, and aircraft reference frame's, 
and finally the' ground-fixed reh'rence frame'. There' is a certain amount of fh'xibilitv in 
delining the vehicle relative to a given reference frame', for instance, the definition of the 
blade shape' in the 1 blade' rc'lc're'iice' frame 1 may be' base'd on zcto sc'ctional twist angle' at the' 
root, the 1 tip, or 75% blade radius. Of course, the rotor e*ollc'clive pitch angle 1 must then be* 
adjust e*d approprial ely. 
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5.2 Harmonic Representations 

Since most of ROTONET is written with the assumption 1 hat the airloads, and thus 
(he noise, is periodic, it is convenient to carry the aerodynamic response information in a 
Fourier series. To relate the Fourier sine/cosine series coefficients to the array form used by 
ROTONET. the following table is given for the translation of analytical Fourier coefficients 
to array elements for an arbitrary quantity A(t): 


liOTO SET form 

Fourier srm.s 

A( 1) 

ft 0 

A( 2) 

ft l 

A (A) 

ft 2 

U.V/2) 

f*\/2-\ 

-U-V/2 + 1) 

ft-S/2 

U.V/2 + 2) 

ft -X/2+1 

M -V - 2) 

ft 

MV- l) 

0-2 

/MA) 

ft - 1 


Table 5.1: ROTONET Harmonic Array Format 
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A(/)= E (5.1) 

ii— —/V 

is the time history in question. Note the explicit use ol the negative Fourier coefficients, 
though the information is already present in the positive coefficients. These harmonic tables 
occur frequently in the higher harmonic airloads modules such as RIX. RRI). and Rid). 

’I he lol lowing table briefly lists the commonly used modules for ROTONET analyst's, 
the inpul/out put flows, and a short description: 


Module 

Inputs From 
Modules 

Outputs to 
Modules 

Purpose 

ABL 

none 

TIN 

atmospheric turbulence properties 

ABS 

ATM 

SFO 

atmospheric absorption coefficients 

ARM 

AIM 

various 

atmospheric properties in parameter form 

AIM 

none 

ABS, SFO. PRO 

atmospheric properties 

EFF 

LEV 

none 

EPNL for each ground observer 

CKO 

SFO 

PRO. PRT 

source/observer geometrv 

IBA 

IBS 

I HI 

inviscid aero force's and moments 

IBL 

IBS. IBA 


viscous force's 

IBS 

none 

IBA, IBL 

reads and interpolates blade shape' 

LEY 

PRO. PRT 

EFF 

frequenev-integrated noise' met t ie s 

LRN 

IBS, IBA. LRP 

MSN. PRT 

compact-chorelwise rotor lone' noise 

LRP 

IBS. IBA, IBL 

LRN. RBN, RLI) 

rolor trim and aerodynamics 

MSN 

LRN. RTN. RBN 

PRT. PRO 

combine's sou ire noise for two rotors 

PRO 

RBN. TIN, MSN 

LEV 

propagates rotor broadband noise' 

PRT 

LRN, RTN. MSN 

LEV- 

propagates rotor tone noise' 

RBN 

LRP 

PRO, MSN 

rotor broadband noise' 

RIN 

Rid). RVVC 

Rid) 

wake- induced velocity 

Rid) 

RIN. RRI) 

RIN. RRI), RTN 

rotor airloads 

RRI) 

Rid). RRI) 

Rid). RRI) 

rigid-blade flapping response 

RTN 

Rid). RRI) 

MSN, PRT 

lull-surla.ee rotor tone noise 

RVVC 

none 

RIN 

rigid wake' geomet rv 

SFO 

AIM 

CEO. PRO. PRT 

steady ic've'l flight trajectory 

STL 

none 

TIN 

streamline distortion ten, sen* 

TIN 

ABL. STL 

PRO 

turbulence ingestion noise' 


Table 5.2: ROTONET Modules 

note that the IBS/IBA/IBL results art' also used by LRE. LRN, R BN. RI.V RRI). Rid), 
and RTN. The following figure illustrates the basic flow of information between modules: 
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Chapter 6 


INPUTS 


6.1 Parameters 

Parameters are variables ol integer, real, logical, or string content, and can either be 
single element or arrays. Their valne(s) can be changed with the PARAM or EVALUATE cont rol 
statements. PARAM statements assign the value of a constant into the parameter. EVALUATE 
statements assign the evaluated value ol an algebraic expression to the parameter. Param- 
eters exist in run-time only; i.e. they can not be stored in a database for use in a later 
RO I ON hi run. Some parameters are changed by Junctional modules: lor example, the 
parameter AO. collective pitch, is set as an initial guess going into module LRP, and contains 
the final predicted value after LHP finishes execution. Other parameters, such as IOUT and 
IPRINT. are not changed by funct ional module's. Examples of syntax can be found in Chapter 
1 and Appendices (’. 1). and E. and detailed discussion of syntax is found in the “Aircraft 
Noise* Prediction Program User's Manual” [2]. 

Note alse> that the' ANOPP Control Language has a rename feature. This allows vou 
te> use 1 a custom parameter name' which is eepiivalenced te> the default parameter name (as 
found in the module's documentation). A typieal situation for this feature' is for a helicopter 
with a main and tail rotor configuration. Since separate LHP runs for each rotor will create 
name conflicts lor parameters such as AO (rotor collective pitch angle in radians), the rename 



ROTONET Primer 


6-2 


feat ure allows new parameter names such as MRAO and TRAO to represent- main rotor and tail 
rotor collective. 

Appendix A contains a global list of default names for parameters commonly used in 
ROTONET analyses. This list expedites development of new input decks, and speeds trace- 
hack of a parameter's value history. During modification or creation of an input deck, the 
proposed name of a new parameter should he checked against the global del an It namelist to 
avoid inadvertent choice of a parameter name already in use. Also, if a parameter value is 
different than that defined by the user, the list allows determination of which modules may 
have overwritt en its value. 

6.1.1 Notes on certain parameters 

The initial guess for rotor collective pitch angle. AO. is used by the LRP module in its 
iterat ive determination of simplified rotor trim. II AO is too far away from the trimmed value. 
LRP will not converge and an error message will be printed to the standard output file. The 
IRL module allows use of boundary layer analysis using flat-plate or curved-wall methods. 
If the drag force table generated by IBL is unrealistic, set the FLAT input parameter to 
.TRUE, and rerun the module. The GLIM parameter used by LRN determines the relative 
error tolerance of the retarded time equation iteration; when in doubt, set it very small to 
avoid inducing phase errors in the predicted time series. In general, the TETCAL and TVFCAL 
models in RBN should be used for typical rotor problems; the- trading-edge blunt ness noise 
and laminar vortex shedding noise models usually contribute little 1o overall noise levels. 

6.2 Members 

f : nit members are data arrays which are fixed and not interpolated, observer locations 
for example. Unit members are created with the UPDATE command. Tables, also called t able 
members, are data base members which are intended to be interpolated, blade shape lor 
example. Their format is known as 1 ype 1 in the ANOPP literature. When a table is 
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creat-exl, the user specifies the kind of interpolation to he used. Tables are created with the 
TABLE cont rol statement. Examples of syntax can he found in Chapter 4 and Appendices (\ 
1), and E. 

Note also that the ANOPP Control Language has a rename feature. This allows you 
to use a custom unit, table, or member name which is equivalenced to the default member 
name (as found in the module's documentation). A typical situation for this feature is for a 
helicopter with a main and tail rotor configuration. Since separate LHP runs for each rotor 
will create name conflicts for members such as IBS (SHAPE), the rename lea t lire allows new 
dat a unit names such as MRIBS (SHAPE) and TRIBS (SHAPE) . These t wo members could also 
be renamed IBS(MRSHAPE) and IBS (TRSHAPE) ; however, since there are many more member 
names than unit and table names, the former scheme requires fewer renames and thus is 
easier and less prone to errors. 

Appendix B contains a global list of default names for members commonly used in KO- 
TONET an ah ses. This list expedites development, of new input docks, and speeds traceback 
of a member's run-time history. During modification or creation of an input deck, the pro- 
posed run-time name of a member should be checked against the global default name list 
to avoid inadvertent choice of a member name already in use. Also, if a member state* is 
different than that defined by the user, the* list allows determinat ion of which module's may 
have* overwritten its value's. 

6.2.1 Notes on certain members 

The SFIELD(PHI) mnnber, use'd to input the lateral distribution of obse'ners on the 
sound bubble', has a unique feature*. If the flight trajectory is parallel to the* A axis in the* 
ground coordinate's, the ground observe*!* array is rectangular with lenigitudinal (along the* 
flight path) rows e>f observers parallel to the A axis, and the* flight path is leved over a flat- 
ground planey them there are 1 constant value's of PHI which will corre'spond to e'mission angle's 
of acoustic "rays" which will re'acli the* longitudinal obse'rver rows precisely. By calculating 
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these emission angles, the number of PHI values can be limited to the number of longitudinal 
observer rows, thus eliminating the need to interpolate the source noise on the bubble in the 
lateral (perpendicular to the flight path) direction. By minimizing the number of observers 
on the sound bubble, the run times will be reduced and accuracy increased for both the 
source noise deck and the propagation deck. 
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EXAMPLE CASE 


] his chapter presents three input decks which, when executed in sequence, predict the 
EPNL tor a level flyover ol a McDonnell Douglas MD-500K helicopter, at 95 knots airspeed. 
3000 pounds gross weight, and 250 feet altitude. Observe that the sequence of deck execution 
is such that the blade aerodynamic properties, rotor trim M airloads calculations, and source 
noise are all calculated before the third deck determines the atmospheric properties as a 
function ol altitude. This is allowable because the flyover altitude is only 250 feet, and 
the ground plane is assumed to be at sea level. I ll is is an acceptable procedure for most 
KOTONHT calculations. How ever, lor flyovers with high altitude above ground level and/or 
Altitude above sea- level, the decks would have to be reconfigured to compute air properties 
(such as density, sound speed, and kinematic viscosity) at the flight altitude of the aircraft. 

1 he process of building input decks can be done in several ways. The beginning 
HO J'ONlvI user is probably best off using existing decks as "templates" for the new decks. 
The danger in t his approach is that re-using various input parameters, unit members, and 
table members blindly can result in incorrect predictions; some thought should be given to 
all inputs. 1 he code' has defaults lor almost all input parameters. For many parameters, 
t hese can be relied on to give good results, but there are many parameters t hat really should 
be customized 1 o t he part icular input deck. The besl source of informal ion for building new 
input decks are t he help files (also called "pseudocode" ). It is recommended t hat when a 
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new version of ROTONET is recievod. the documentation for all modules pertinent to RO- 
TONET be extracted and printed out . That printout is t he most frequently used document 
when advanced users run ROTONET. 

7.1 Blade Properties Deck 

The input deck listed in Appendix (' executes the IBS. IBA. and IBL modules, creating 
a library file which contains t he blade geomet ry and aerodynamic characteristics for t he main 
rotor. 

The IBA (LIFT). IBA (MOMENT) . and IBL (DRAG) tables should be plotted out to chock for 
“fiver" points and to make sure the lift and drag curve slopes are reasonable. Most airfoils 
used on helicopter rotors behave like traditional NACA sections. For instance, the lift curve 
slope should be somewhat less than 2tt. A good place to check the curves is in Abbott iV von 
Doenhoff [27]. Note that the ROTONET blade properties module's do not have sophist icated 
stall modeling, so comparisons beyond t he actual stall angle' will not make sense'. A stall 
cutoff parameter is available to limit the ROTONk/1 -predict f'd lift to a known stall value. 
The interpolation capability's of the blade shape modules ate limited. Thus, il the user 
wants a hue grid at the leading edge', plenty of airfoil leading t'dge points should be used in 
GE0M( BLADE). While increasing resolution also increases run time lor black' properties decks, 
it is generally wise to have fine resolution, since these decks arc' only run once per rotor. 
One can always use a coarser grid in the actual noise calculations. I he (iRII)(XI2) member 
represents t he chordwise angular coordinate in the Joukowski transform plane. The values 
run from 0.0 at the upper trailing edge, to 0.5 at t he leading edge, to 1.0 at the lower trailing 
edge. The grid shown in this demonstration deck is one developed from experience for full- 
surface modeling wit h RTN. It is recommended as a, st andard grid for most calculat ions. Note 
that the spanwise distribution used in GEOM(BLADE) will be the one output into GRID(XIl). 
If the blade' has unusual feature's such as nonlinear twist , rapidly changing airfoils, breaks in 
plan form etc., then the spanwise st at ions used in GEOM (BLADE) must be enough to adequately 
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describe the blade. Note that if the GRID(XIl) member desired for the noise computations 
is not the same as that produced by IBS (which is quite likely), the user should rename 
GRID(XIl) on the "EXECUTE IBS"" command or put the command PURGE GRID $ alter the 
"EXECUTE IBS statement. Also, if the user gets a "zero determinant" error from IBS or a 
"laminar separation" message from IBL, the solution is to get a liner blade shape' definition 
in GEOM(BLADE) and a finer GRID(XI2) distribution. 

7.2 Source Modeling 

The deck listed in Appendix 1) calculates the main rotor aerodynamic performance', tone 
noise, and blade' broadband self-noise, using the modules LHP. LRN, and RBN. 

Like 1 any helicopter trim program, t here are' certain t hings to look for in a rotor airloads 
prediction, for the given flight condition, the gross behavior like collective pitch and tip- 
path-plane' angle' should be compared with other trim codes or experimental data. The cyclic 
pitch inputs and flapping response is the next level to check, again with comparison to other 
predictions or empirical data. Finally, the Mach number and angle-of-attack distributions 
should be looked at to ensure there' are' no nonsensical aerodynamic conditions being output 
(eg. a overly large' angle-of-attack at the blade tip in forward flight). The GRID(XIl) 
member usee l in the' above 1 ele'ek has been devebped from experience with various aircraft 
ease's. Note' that, if the eolle'ctive pitch and/or angle-of-attack iterations fail to converge, 
putting a realistic initial guess in for the parameter AO tenuis to help. 

The predicted tone noise* at the sound bubble can be' checke'd by plotting the pressure 
time' histories and narrowband spectra, and comparing them with wind tunnel data for simi- 
lar flight conditions and observer location. Several publications are available' with measured 
data [28. 29. ♦?()]. (lenerally speaking, the* thickness noise component should be maximum 
in the rotor plane, and is highly directional for a conventional rotor. Loading noise should 
be maximum below the* rotor in forward flight and is le'ss directional (excepting BVI). It has 
been lounel that NTIME should be' four time's greater than NHARM. Ofteui. a factor of two is 
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used due to the Nyquist criteria, hut in practice, better results are obtained with a factor of 
four. 

7.3 Flight Trajectory & Propagation 

The input deck listed in Appendix E tabulates atmospheric properties, aircra.lt flight 
pat h, and source/observer geometries, and propagates the source noise at the “sound bubble' 
to the ground observers. Modules A US, SKO. CEO, PRO, PEI . LEV, and EEE are used. 
The end result of this deck is time histories ol metrics such as PNLT. as well as the EPNL 
for the entire flyover. 

The user should check the results from SEO to ensure that the time step, initial, and final 
flight times will result in a flight history which brackets the desired inilia] and filial observer 
elevation angles. If one relies on the defaults, the flyover may not include flight times which 
contribute substantially to EPNL. The easiest wav to check this is to look at the elevation 
angle printout for each observer. For rotorcraft , this angle* should begin at roughly 10 or 20 
degrees, and end likewise. The 4-octave band center frequencies in SFIELD(FREQ) should 
be consistent with those used downstream in the propagation and integration modules. For 
purposes of computing PNLT and EPNL, the 24 standard bands (with center frequencies 
starting at 50 Hz and ending at 10 KHz) should be used. However, to correctly account for 
the low frequency content of main rotor noise, the set of bands should be extended down to 
the 12.5 Hz center frequency. 

The best bet for checking the predicted noise propagated to ground observers is to 
compare the predicted spectra with empirical data. Unfortunately, the database* of flyover 
measurements is limited to a lew vehicles [51, 42, 44. 44]. However, these* tests do allow a 
certain amount of trend information to be checked. The time history of PNLT should lx* 
checked to ensure that the start and end times of tlu* flyover are large* enough that the start 
and end PNLT values are more than 10 dB lower than the peak PNLT. This is necessary 
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to correctly coni])utc hPNL. Also, ensure that 1 1 k' SFIELD(FREQ) array used is identical to 
that used by the (!h() module as well as the LRN and RUN module's. 
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Global Parameter List 


Parameter 

Form 

Created By 

AO 

re'al 

user. LRP 

AOL 

real 

ERE. RLI) 

a or 

re'al 

RED 

Al 

re'al 

ERR 

A 1 L 

re'al 

ERR. REI) 

All’ 

real 

RLI) 

A A 1 

re a al 

user 

A A 2 

real 

user 

A MSI NT’ 

integer 

user 

ABSORB 

logical 

user 

A ERROR 

re'al 

user 

AKL 

real 

ARM 

AERDR 

re'al 

ERR 

A LRU A 



aertrr 

re'al 

RED 

AMI 

real 

ARM 

AXCEES 

re'al array 

user 

AX! 

re'al 

ARM 

arrend 

logical 

user 

AW 

revd 

user 

Bl 

re'al 

LHP 

B1L 

real 

LHP, RLD 

mr 

re'al 

HL1) 

BBB 

re'al 

user 

BETA 

real array 

user 

BN 

integer 

user 

BT 'IP 

logical 

user 


Description 

collective pilch angle 
coiling angle 

updated collective pitch angle 
lateral cyclic pitch angle 
longitudinal cyclic flap angle' 
updated lateral cyclic pitch angle 
rotor pitch angle relative to body 
rotor pitch angle relative to body 
number of int egral ion steps in atmospheric 
absorpt ion calculation 

acti vates/deact ivat.es atmospheric absorption 

error tolerance lor streamline integration 

thermal conductivity 

hub plane' angle' 

aircraft angle'-of-at t ack 

rote>r tip-pa.th-pla.ne 1 angle 1 

dynamic vise’osity 

array of direct ivity angle's 

kinemat ic viscosity 

flag for SI () to create 1 a nenv flight profile 1 or 
append to the existing one in FLI(PATII) 
reference 1 area e>f aircraft 

longitudinal cyclic pitch 
lateral cyclic flap angle 
undated longitudinalcyclic pitch 
blaele 1 tip loss factor 

0-th, 1st, 2nd rigid flapping harmonics (can 
be assembled from AOL. AIL, H1L. see 
exam])le eh'cks) 
number of blades 

flag to choose blunt /not blunt blade 1 tip 
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Parameter 

Form 

Created By 

Description 

r 

real 

user 

blade chord 

GA 

real 

user. ARM 

sound speed 

('A HAH 

real 

ARM 

average sound speed 

Cl. MAX 

real 

user 

maximum lift coefficient allowed 

(■OH 

real 

user 

incoherence coefficient 

GO 

real 

LRR. RLD 

torque coefficient 

CT 

real 

LRR. RLD 

thrust coefficient. 

CTK 

real 

user 

characteristic t ime const ant 

ELDB 

real 

user 

dB threshold below max PNFT for 
EPNL calculation 

DKLH 

real 

user 

altitude increment 

DELF 

real 

user 

frequency bandwidt h 

DELMACII 

real 

user 

Mach number limit 

DELTA 

real 

user 

flap setting (irrelevant to ROTONET) 

DELHI 

real 

user 

maximum polar directivity angle limit 

DIRECT 

logical 

user 

flag to use input ./computed reception times 

1) ROTOR 

real 

LRR 

rotor drag force 

DS 

real 

user 

distance 1 between streamlines 

DTAIL 

real 

user 

distance between main k tail rotor hubs 

Dll ME 

real 

user 

ground reception time increment 

e 

real 

user 

rotor blade flap hinge 1 offset 

flag to choose EPNL or SPL output 

EFFORT 

integer 

user 

ENGNAM 

Hollerith 

user 

sets XXX pari of FLl(FLIXXX) member name 

EMF 

logical 

user 

flag to use of external surface* pressures 

ERSLON 

real 

user 

convergence* limit lor IBA 

F 

real 

user 

vehicle* flat plate drag area 
larfield distance* for FIN calculations 

FF 

real 

user 

FI LX AM E 

Hollerith 

user 

name of file for contour data 

Ell 

real 

user 

blade moment of inertia 

ELAT 

logical 

user 

flag to choose* boundary layer model 

ETZERO 

real 

user 

fraction of eventual tip vortew streuigth 
initially rolled up into tip vortex 

C 

real 

user 

geostrophic wind spevd 

GAMMA 

real 

user 

specific heat ratio 

OCT 

real array 

user. RKI) 

constants used in RR1) calculations 

GEO ERR 

integer 

CKO 

flag dened.ing success of (JEO calculations 

GLIM 

real 

user 

erre>r toleranee* for regarded time calculation 

GROUND 

logical 

user 

flag to enabh'/disable* ground e*ffects 

ill 

real 

user 

ground level alt it ude 

HE 

integer 

integer 

user 

counter of first harmonic in TIN analysis 

HL 

user 

counter of last harmonic in TIN analysis 

USER 

real 

user 

shape factor at se*paration point 

HUMID 

real 

A PM 

humidity (se*e* A PM ) 
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Parameter 

Form 

Created By 

Description 

I AERO 

integer 

user 

LRP/RLl) airloads selector for HHN 

I AWT 

logical 

user 

enable/ disable dHA output 

IBLFRCT 

Hollerith 

user 

member name prefix for IBL(YYYNNN) 

I RAPHES 

Hollerith 

user 

member name prefix for IBA(XXXNNN) 

I('OORI) 

integer 

user 

source/observer coordinate selector 

IDWI 

logical 

user 

enable/disable dBI) output 

III 

integer 

user 

turbulence selector 

I MODE 

integer 

user 

MSN mode selector flag 

IMPROY 

logical 

user 

blade properties format selector 

IOP 

integer 

user 

flight condition/observer format selector (RBN) 

10 IT 

integer 

user 

contour data selector ((’NT) 
coordinate system option flag (MSN) 

IOSPL 

logical 

user 

enable/disable OASPL output 

10 IT 

integer 

user 

printout format selector (LEV. PRO. PRT. HHN) 
output member format selector (LRN. H I N) 

IOETPUT 

integer 

user 

output option flag 

IP XL 

logical 

user 

euiable/disable 1 PNL out put 

IPX LI 

logical 

user 

enable 1 / disable PNL'T output 

1 PHILS 

ini eger 

user 

plot file 1 out put- selector 

IPRINT 

integer 

user 

input /out put print selector 

IHOLL 

integer 

user 

tip vortex rollup extent- index 

is\ i: 

i tit eger 

user 

atmospheric absorption selector 

is r h 

integer 

user 

enable/disable streamline file output 

ISVM 

logical 

user 

symmetric ground observer array flag 

MAIL 

integer 

user 

main/tail rotor trim analysis selector 

ITEM E 

ini eger 

user 

enable/disable deformation tensor file 1 output 

UMAX 

integer 

user 

iteration limit for retarded time ('((nation 

HAILS 

I Inherit h 

user 

unit s select or 

J 

integer 

SFO 

initial flight time step counter 

L 

real 

user 

Monin-Obukhov stability length 

LBLCAL 

logical 

user 

enables/disables laminar boundary layer noise* model 

LOSLOP 

rc'al 

user 

lift curve slope 1 

MACHO 

real 

user 

turbulence convective Mach number 

MACH HE 

real 

user 

rotor hover tip Mach number 

MASSAC 

real 

user 

reference mass of the aircraft 

MEMSEM 

Hollerith 

user 

unit iV member names of broadband noise 1 
tables to be summed 

MEM SI AIN 

Hollerith 

user 

unit V member name's of narrowband noise 1 
tables to be summed 

MIS 

real array 

user. HR I) 

flapwise moment array 

ME 

real 

user 

rotor advance 1 ratio 

MVVB 

real 

user. LHP 

blade flapping moment of inertia 

MZ 

real 

user 

aircraft, flight Mach number 
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Parameter 

Form 

Created By 

Description 

N AWT 

logical 

user 

enable/disable A- weighting of narrowband data 

NBANI) 

integer 

user 

nuber of 1/3-octave band sub-bands 

N BLADE 

integer 

user 

number of blades on rotor 

NCOFFF 

integer 

user 

number of element s in BETA and ZETA arrays 

NDWT 

logical 

user 

enable/disable D- weighting of narrowband data 

NFRR 

logical 

many modules 

denotes error status of job stream 

N II ARM 

integer 

user 

number of harmonics to output; 
acoustic pressure (LRN, RTN, MSN), 
azimuthal (KIN. RL1). RRI)) 

NIK) 

integer 

user 

number of altitudes at which to 
output atmospheric quantities 

NISO 

integer 

user 

number ol streamlines for which 
trajectory is to be computed 

NORUFR 

integer 

user 

order of sine series interpolation of blade shape 

NOS PL 

logical 

user 

enable/disable calculation of OASPL 
from narrowband noise* data 

NPSI 

integer 

user 

number of azimuth angles for LRP. RW(I 

NR 

integer 

user 

number of radial stat ions for TIN noise* 

NSS 

integer 

user 

max time steps used in streamline* integration 

NREV 

integer 

user 

number of wake recs used for wake inflow 

NT 

integer 

user 

number of azimuth stops for TIN noise* 

NTIMF 

integer 

user 

number of time steps used in acoustic 
analysis (LRN, RTN) 

NFM 

integer 

user 

number of frequencies in TIN noise* 

NX 

integer 

user 

number of X points reemireel for 
('artesian moele calculalieui e)f streamliners 

NY 

integer 

user 

number of Y points 

NZ 

integer 

user 

number of X point s 

OMFCJA 

real 

user 

re)tor rotat ional spevd 

OMR 

real 

user 

rotor tip speeel 

PI 

real 

user 

atmospheric pressure* at ground leyel 

PARNAM 

Hollerith 

user 

array e>l parameter names te> be output- 

pit; 

Hollerith 

user. SIX) 

initial landing gear position 
file name for LRN ple>t file* 

PLT1I) 

Hollerith 

user 

PRFS 

real 

A PM 

atnmspheric pressure* (se*e* ARM) 

P RON PM 

integer 

user 

noise* table number (see LRN(XXXNMN)) 

PROSUM 

Hollerith 

user 

names e>f broadband noise* members to 
be summe'd be*fore propagation 

PROTIMF 

Hollerith 

user 

member name* suffix of FLI(FLIXXX). 
alse> use*el as member name prefix of PRO output 

prttimf 

Hollerith 

user 

member name* suffix of FLI (FLIXXX) . 
also used as me*mbe*r name* j)iefix e>f PR I output 

PSIO 

real 

user 

initial rote>r azimuth angle* 
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Parameter 

Form 

Created By 

H 

re'al 

user 

HO 

real 

user 

HI 

re'al 

user 

H2 

real 

user 

HBAPRKS 

Hollerith 

user 

HKC'POL 

Hollerith 

user 

HKHHOH 

re'al 

user 

HI 10 A 

re'al 

user, ARM 

HHOACA 

re'al 

ARM 

H INI 

re'al 

user 

HL 

re'al 

user 

HNDTIP 

logical 

user 

HOTLFFT 

logical 

user 

HS 

re'al 

LRN 

HX 

re'al 

user 

SCHNNN 

integer 

user 

SCHXXX 

Hollerith 

user 

SIDLIN' 

logical 

user 

SIGMA 

re'al 

user 

SOL 

nvd 

user 

SHUCK 1 

real 

user 

SHHOF 

real 

user 

SIMM 

real 

user 

STATIC 

logical 

user 

STOP 

re'al 

user 

SUHFACK 

Hollerith 

user 

I’A 

re'al 

ARM 

I AIM 1 1 

re*al 

LRR 

TB2CAL 

logical 

user 

TBLCAL 

logical 

user 

TFBCAL 

logical 

user 

IKTOAI. 

logical 

user 

II 

re'al 

SFO 

ruin \ 

re'al 

us<'r 

ri i ho r 

real 

user 

rnvv 

real 

user 

n 

real 

user 

I LO 

re'al 

user 

TSTFP 

re'al 

user 

TV FOAL 

logical 

user 

VO 

revd 

US('r 

VF 

re'al 

user. SFO 

VI 

re'al 

user, SIX) 

VI XF 

re'al 

user 


Description 

blade radius 

source/observer arc radius (not used) 
rotor 1 radius 
rotor 2 radius 
not used 

polar/C artesian streamline coordinate flag 

error tolerance in STL streamline integration 

ambient air density at aircraft 

characteristic impedence 

Reynold's number ( = R*(A/ANE) 

rotor lag hinge offset 

round/s(|ua.re tip flag 

left/right -hand rotation flag 

sound bubble radius, dimensional 

sound bubble radius, nondimensional 

member name suffix for XXXNNN 
member name prefix for XXXNNN 
flag to choose between observer array 
on ground plain' or arc 
specific flow resistance of t he ground 
rotor solidity rat io 
| unburst wake con' radius 
burst wake con' radius ratio 
initial flight time 
fixed/moving noise source' flag 
final flight t ime 
hard/soft ground surface flag 

ambient temperature' (see ARM) 
required tail rotor thrust 
enable/disable Brooks' TEB-VS noise 1 model 
enable/ disable' Schlinker/ Amiet 
TBL-TE noise 1 model 

enable/disable Grosveld's TEB-V S noise 1 mode*! 
enable'/disable Brooks' TBL-TE noise mode 1 ! 
final flight time 

geographic latitude (note 1 name' conflict 
with SFIELD (THETA) member) 
throttle 1 setting (not relevant for ROTONET) 
flight climb angle 1 
initial flight time' 
landing gear reset time 
time step for flight trajectory 
enable'/disable' Brooks' t ip vortex 
formation noise model 

average' wake-induced inflow velocity 
final flight velocity 
initial flight velocity 
vertical climb velocity 
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Parameter 

Form 

Created By 

Description 

W 

real 

user 

aircraft weight, coefficient- (RLD) 

WB 

real 

user 

total blade weight 

WEIGH T 

logical 

user 

weighting opt ion for GEOM (BLADE) 

WKPRNT 

logical 

user 

wake distortion/position print flag 

WT 

real 

user 

aircraft weight coefficient (LRP) 

XI 

real 

user 

minimum X for streamline calculat ions 

X2 

real 

user 

maximum X for streamline calculations 

X ARRAY 

real array 

user 

longitudinal ground observer distribution 

XF 

real 

user. SFO 

final longitudinal flight position 

XI 

real 

user 

initial longitudinal flight position 

XMAX 

real 

user 

streamline domain limit (see STL) 

XMIN 

real 

user 

streamline domain limit (see STL) 

Y1 

real 

user 

minimum Y for streamline calculations 

Y'2 

real 

user 

maximum Y lor streamline calculations 

YARRAY 

real array 

user 

lateral ground observer distribution 

YF 

real 

user. SFO 

final lateral flight position 

YI 

real 

user 

initial lateral night position 

YMAX 

real 

user 

streamline domain limit (see STL) 

YMIN 

real 

user 

streamline domain limit (see STL) 

Z 

real 

user 

height above ground 

zo 

real 

user 

ground surface roughness height 

Z1 

real 

user 

minimum Z for streamline calculations 

Z2 

real 

user 

maximum Z for streamline calculations 

ZETA 

real array 

user 

rigid blade lead/lag angle harmonics 

ZF 

real 

user, SFO 

final vertical flight position 

ZGR 

real 

user 

altitude of ground plane above sea level 

ZI 

real 

user 

initial vertical flight position 

ZOPT 

integer 

user 

flight trajectory flag 

zs 

rca 1 

user 

source altitude (for nonmoving source 1 ) 

ZSLOPE 

logical 

user 

flag for blade shape interpolation 

ZM AX 

real 

user 

streamline domain limit (see STL) 

ZMIN 

real 

user 

streamline domain limit (see STL) 
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APPENDIX B 


Global Member List 


Member Form Created By Description 

ATM (A AC) table member ABS atmospheric absorption coefficients 

ATM(IN) unit member user temperature A- .humidity table 

A I Mi l MOD) table member APM. ATM air properties table 

also avaiable from LIBRARY AVRK 
as ATMOS(STHD) 

BFAM(STFFR) table member user arra.v associated with beam 

steering analysis (not used) 

BLAI)K(TRIP) unit, member user boundary layer flip location array 

I'd' 1‘ ( FPNL) unit member KFF FPNL at ground observe*!' arrav 

FFF(SFF) unit member FI T’ SFL at ground observer arrav 

FIS(VFb()(TrV) table member user rotor inflow from external source 

FMP(PRFS) table* member user blade* surface* pressure coefficient 

FMP(SHFAH) table* member user blade* surlae’e* shew e‘oe*f 1 i cie*i it 

FNV(COFF) table member use*r wake* geometry e*nvele)pe functions. 

available from' RYYCLIB.WRK 

FLI(FLIXXX) unit member SFO senure* t ime*/(light .condition table* 

FLI(PATIl) unit member SFO flight trajectory table* 

CFOfCFOM ) unit member CFO CFOM = BODY. WIND. XXXX 

OFOM(BLAI)F) unit member user prop/ rotor blade geomet ry 

CRI)(( 'OORD) unit, member CRI) ground observe*!* array 

CRID(PSI) unit member user table of e.*s for FRP airloads 

CRID(XIl) unit member user blaele* radial grid points 

CRII)(XI2) unit member user blade cliordwise grid points 
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Parameter 

Form 

Created By 

IBA(AEROCTR) 

unit mei)ib(T 

IB A 

IBA( ALPHA) 

unit member 

user 

IBA(INCOMP) 

unit member 

IB A 

IBA(LIFT) 

unit member 

IBA 

IBA(MACIl) 

unit meml)er 

user 

I BA( MOMENT) 

unit member 

IBA 

IBA(STACi) 

unit member 

IBA 

IBA(VYVNNN) 

unit member 

IBA 

IBA(ZERO) 

unit member 

IBA 

IBL(l)RAG) 

unit member 

IBL 

IBL(THKK) 

unit member 

IBL 

IBL(YYYNXN) 

unit member 

IBL 

IBS(El)GTIIK) 

unit member 

IBS 

IBS( SHAPE) 

unit member 

IBS 

IBS(SLPCRD) 

unit member 

IBS 

IBS(SLPSPN) 

unit member 

IBS 

IBS(SPAN) 

unit member 

IBS 

LEV(AWGT) 

unit member 

LEV 

LEV(PNLT) 

unit member 

LEV 

LRP( FORCES) 

unit member 

LRP 

LHP( PER FORM) 

unit member 

LHP 

LRP(STRNTII) 

unit member 

LRP 

OBSERV(COORD) 

unit member 

user. 

CRD. 

MSN 

RIN(QMAT) 

unit member 

R 1 X 

RIN( VELOCITY) 

unit member 

RIN 

R LI )( FORCES) 

unit member 

RL1) 

RLD( GAMMA) 

unit member 

RLI) 

RLI)(M OTION) 

unit member 

RLD 

RRI)( ANGLES) 

unit member 

RRD 

RRD(MASS) 

t able member 

user, 

LRP 

RSD(BEND) 

table member 

user 

REA(LIET) 

t able member 

user 

RWC(WKI)IST) 

unit member 

RWG 





Description 

aerodynamic (enter distribution 
angles-of-attack for 1BA calculation 
incompressible pressure coefficients 
section lift coefficients 
angles-of-attack for 111 A calculation 
section moment coefficients 
stagnation point distribution 
blade surface pressure 
zero lift angle-of-attaek dist il 
section drag coefficients 
boundary layer thickness 
skin friction table 
trailing edge t hickness dist il 
blade shape t able 

chordwise derivative of blade shape 
spanwise derivative 1 of blade shape 
blade cross-sectional characteristics 

(1BA time history at ground observers 
PNLT time history at ground observers 
blade airloads tables 
blade aerodynamic angles A’ velocities 
max bound circulation table 


see GK1)( COORD) 
see MSMOBSK OBS2) 
ground observer locations (CEO), 
nubble observer locations in 
Cartesian form 
(LBN. RTN. RUN. TIN) 

wake influence coefficients 
wake inflow at rotor disk 
blade airloads table 
blade bound circulation table 
section Mach number A r 
angle-of-att ack tables 
rigid flapping harmonic table 


blade mass dist ribution 
blade bending harmonics 
unsteady aero corrections 
wake geometrv distortion table* 
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Parameter 

Form 

Created By 

Description 

SIIP(('OS(’) 

table 1 member 

user 

from RWCLIB.VVRK 

SIIP(SINC) 

table 1 member 

user 

from RWCLIB.VVRK 

SlIPfZKROC) 

table 1 member 

user 

from RWCLIB.VVRK 

SFIKLD(KRKQ) 

unit member 

user 

I/d oc tave 1 band renter frequencies 

SFIKLD(PIII) 

unit, member 

user 

sound bubble observer azimuth angles 

SFIKLD(TIIKTA) 

unit member 

user 

sound bubble 1 observe-M* polar angle's 

STLTIX( ABLOTl ) 

unit member 

ABL 

at mosplmrie* turbulence 1 propeut ie's 

stltin(stloti) 

unit member 

STL 

rotor inflow distortion tensor 

TIM(HISTl) 

unit member 

LRN. 

see LRN(TIMK) 



RTN 

see RTN ( TIM K) 

TI.\1(IIIST2) 

unit member 

LRN, 

see LRN(TIMH) 



RTN 

see RTN(TIMK) 

TIP( AXCFF) 

table member 

user 

tip angle' correct ion factor 

XXXXXX(YYYYYY) 

unit member 

PRO 

see PRO(PRKS) 

Y Y Y Y Y Y ( X X X N N N ) 

table 1 member 

LRN. 

see LRN(XXXNNN) 



RTN. 

see RTN(XXXNNN) 



MSN. 

see MSN(XXXXNX) 



RBN, 

see RBN(XXXXNX) 



TIN 

see TIN(XXXNNN) 

ZZZZZZ(WWWWWW) 

unil member 

PRT 

SFH PRT(PRFS) 
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APPENDIX C 


Example Deck 1 


ANOPP JECH0= . TRUE . JL0G=. FALSE. NLPPM=55 $ 
STARTCS $ 


############### ########################################################## 

File ; demoOl . inp 

Author; Devon Prichard 
Level; 03/02/11 

Purpose , This input deck computes the blade geometry and aerodynamic 
characteristics of the MD-500E main rotor. 

Input Libraries; none 
Output Libraries; DEM001.WRK 

Modules used; 

IBS - improved blade shape module 

IBA - improved blade aerodynamics module 

IBL - improved blade boundary layer module 

!! WARNING !! this demo produces a large output file! think 
twice before sending output file to printer... 

########################### ############################################## 

Setup Parameters 

#################### ##################################################### 

Define parameters unique to this rotor, otherwise use defaults; 


PARAM R - 13.205 $ blade radius in ft 
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('-> 


PARAM I UNITS = 7HENGLISH $ use english units 

PARAM IPRINT =3 $ print both input and output data 

$ 

EVALUATE RINF = 0.23296E+08 * R * 0.3048 

$ COMPUTE SCALING REYNOLD'S NUMBER 

$ 

$ 

$ ######################################################################### 

$ Setup Members 

$ ######################################################################### 

$ 

$ Spanwise and chordwise grids defined in the following members; 

$ note that span is defined relative to blade radius, and chord 
$ is defined in the elliptical coordinate system, 

$ where the upper TE is theta=0, LE is theta=0.5, and 
$ lower TE is theta=1.0, relative to 2*PI. 

$ 

$ Note that the 500E blade, with basically rectangular planform, 

$ linear twist, and constant cross-section, does not require 
$ a very detailed blade definition. Blades with tip shapes, 

$ airfoil changes, or nonlinear twist distributions need finer 
$ descriptions. With increasing resolution comes increased 
$ CPU time to execute IBS-IBA-IBL, so use advisedly. 

$ 

$ 

UPDATE NEWU=GRID S0URCE=* $ 

-ADDR 0LDM=* NEWM=XI1 F0RMAT=4H*RS$ $ 

0.200 0.400 0.600 0.700 0.750 0.80 0.850 

0.900 0.925 0.950 0.975 1.000 $ 

-ADDR 0LDM=* NEWM=XI2 F0RMAT=4H*RS$ $ 

0.00 0.025 0.05 0.10 0.15 

0.20 0.25 0.30 0.35 0.40 

0.41 0.42 0.43 0.44 0.45 

0.46 0.47 0.48 0.49 0.50 

0.51 0.52 0.53 0.54 0.55 

0.56 0.57 0.58 0.59 0.60 

0.625 0.65 0.70 0.75 0.80 

0.85 0.90 0.95 0.975 1.00 $ 

END* $ 

$ 

$ 

$ The values of sectional Mach number and angle-of -attack in degrees, 

$ at which the blade section aerodymamics are computed, are given here; 

$ 

$ 

UPDATE NEWU=IBA S0URCE=* $ 

-ADDR 0LDM=* NEWM=MACH F0RMAT=4H*RS$ $ 

0.1 0.3 0.5 0.7 $ 
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-ADDR OLDM=* NEWM=ALPHA F0RMAT=4H*RS$ $ 
-6.0 -3.0 0.0 3.0 6.0 $ 

END* $ 


Define the blade geometry. Section is the NACA 0015. Note that 
the trailing edge tab is not currently modeled. Also recall that, 
unlike the usual rotor notation, the blade twist is defined relative 
to zero pitch at the blade root. 


UPDATE NEWU=GE0M S0URCE=* $ 

-ADDR 0LDM=* NEWM=BLADE F0RMAT=0 $ 
5 $ five spanwise stations 

5 $ all same airfoil section 


0.00 

-0.0106 

0.000 

0.043 

0.025 

0.00 

107 

106 

$ 

0% 

station 

0.25 

-0.0106 

-0.000 

0.043 

0.025 

-2.25 

107 

106 

$ 

25'/, 

station 

0.50 

-0.0106 

-0.001 

0.043 

0.025 

-4.50 

107 

106 

$ 

50"/, 

station 

0.75 

-0.0106 

-0.001 

0.043 

0.025 

-6.75 

107 

106 

$ 

75% 

station 

1.00 

-0.0106 

-0.002 

0.043 

0.025 

-9.00 

107 

106 

$ 

100% 

station 


1.00000 .00158 $ 

.99000 .00332 $ 

.98000 .00504 $ 

.97000 .00674 $ 

.96000 .00842 $ 

.95000 .01008 $ 

.94000 .01172 $ 

.93000 .01334 $ 

.92000 .01495 $ 

.91000 .01653 $ 

.90000 .01810 $ 

.89000 .01964 $ 

.88000 .02117 $ 

.87000 .02269 $ 

.86000 .02418 $ 

.85000 .02566 $ 

.84000 .02712 $ 

.83000 .02856 $ 

.82000 .02999 $ 

.81000 .03140 $ 

.80000 .03279 $ 

.79000 .03417 $ 

.78000 .03552 $ 

.77000 .03687 $ 

.76000 .03819 $ 

.75000 .03950 $ 

.74000 .04080 $ 

.73000 .04207 $ 
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.72000 

. 04333 

$ 

.71000 

.04457 

$ 

.70000 

.04580 

$ 

.69000 

.04701 

$ 

.68000 

.04820 

$ 

.67000 

. 04937 

$ 

.66000 

.05052 

$ 

.65000 

.05166 

$ 

.64000 

.05277 

$ 

.63000 

.05387 

$ 

.62000 

.05495 

$ 

.61000 

.05600 

$ 

.60000 

.05704 

$ 

.59000 

.05806 

$ 

.58000 

.05905 

$ 

.57000 

.06003 

$ 

.56000 

.06098 

$ 

.55000 

.06191 

$ 

.54000 

.06281 

$ 

.53000 

.06369 

$ 

.52000 

.06454 

$ 

.51000 

.06537 

$ 

.50000 

.06618 

$ 

.49000 

.06695 

$ 

.48000 

.06770 

$ 

.47000 

.06841 

$ 

.46000 

.06910 

$ 

.45000 

.06976 

$ 

.44000 

.07038 

$ 

.43000 

.07097 

$ 

.42000 

.07153 

$ 

.41000 

.07205 

$ 

.40000 

.07254 

$ 

.39000 

.07298 

$ 

.38000 

.07339 

$ 

.37000 

.07376 

$ 

.36000 

.07408 

$ 

.35000 

.07436 

$ 

.34000 

.07459 

$ 

.33000 

.07478 

$ 

.32000 

.07491 

$ 

.31000 

.07499 

$ 

.30000 

.07502 

$ 

.29000 

.07499 

$ 

.28000 

.07491 

$ 

.27000 

.07476 

$ 

.26000 

.07455 

$ 

.25000 

.07427 

$ 
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.24000 

.07391 

$ 

.23000 

.07349 

$ 

.22000 

.07298 

$ 

.21000 

.07239 

$ 

.20000 

.07172 

$ 

. 19000 

.07095 

$ 

. 18000 

.07008 

$ 

. 17000 

.06911 

$ 

. 16000 

.06802 

$ 

. 15000 

.06681 

$ 

. 14000 

.06547 

$ 

. 13000 

.06399 

$ 

. 12000 

.06235 

$ 

.11000 

.06054 

$ 

. 10000 

.05853 

$ 

.09000 

.05631 

$ 

.08000 

.05384 

$ 

.07000 

.05108 

$ 

.06000 

.04797 

$ 

.05000 

. 04443 

$ 

.04000 

.04035 

$ 

.03000 

.03550 

$ 

.02000 

.02950 

$ 

.01750 

.02772 

$ 

.01500 

.02580 

$ 

.01250 

.02367 

$ 

.01000 

.02130 

$ 

,00750 

.01856 

$ 

,00500 

.01527 

$ 

00250 

.01090 

$ 

00000 

.00000 

$ 

00250 

-.01090 

$ 

00500 

-.01527 

$ 

00750 

-.01856 

$ 

01000 

-.02130 

$ 

01250 

-.02367 

$ 

01500 

-.02580 

$ 

01750 

-.02772 

$ 

02000 

-.02950 

$ 

03000 

-.03550 

$ 

04000 

-.04035 

$ 

05000 

- . 04443 

$ 

06000 

-.04797 

$ 

07000 

-.05108 

$ 

08000 

-.05384 

$ 

09000 

-.05631 

$ 

10000 

-.05853 

$ 

11000 

- .06054 

$ 
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12000 

- .06235 

$ 

13000 

- .06399 

$ 

14000 

-.06547 

$ 

15000 

-.06681 

$ 

16000 

-.06802 

$ 

17000 

-.06911 

$ 

18000 

-.07008 

$ 

19000 

-.07095 

$ 

20000 

-.07172 

$ 

21000 

-.07239 

$ 

22000 

-.07298 

$ 

23000 

-.07349 

$ 

24000 

-.07391 

$ 

25000 

-.07427 

$ 

26000 

-.07455 

$ 

27000 

-.07476 

$ 

28000 

- .07491 

$ 

29000 

-.07499 

$ 

30000 

-.07502 

$ 

31000 

-.07499 

$ 

32000 

-.07491 

$ 

33000 

-.07478 

$ 

34000 

-.07459 

$ 

35000 

-.07436 

$ 

36000 

-.07408 

$ 

37000 

-.07376 

$ 

38000 

-.07339 

$ 

39000 

- .07298 

$ 

40000 

-.07254 

$ 

41000 

- .07205 

$ 

42000 

-.07153 

$ 

43000 

- .07097 

$ 

44000 

- .07038 

$ 

45000 

- .06976 

$ 

46000 

- .06910 

$ 

47000 

-.06841 

$ 

48000 

-.06770 

$ 

49000 

- .06695 

$ 

50000 

-.06618 

$ 

51000 

-.06537 

$ 

52000 

-.06454 

$ 

53000 

-.06369 

$ 

54000 

-.06281 

$ 

55000 

- .06191 

$ 

56000 

-.06098 

$ 

57000 

- .06003 

$ 

58000 

- .05905 

$ 

59000 

-.05806 

$ 
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.60000 

-.05704 

$ 

.61000 

-.05600 

$ 

.62000 

-.05495 

$ 

.63000 

-.05387 

$ 

.64000 

-.05277 

$ 

.65000 

-.05166 

$ 

.66000 

-.05052 

$ 

.67000 

-.04937 

$ 

.68000 

-.04820 

$ 

.69000 

-.04701 

$ 

.70000 

-.04580 

$ 

.71000 

-.04457 

$ 

.72000 

-.04333 

$ 

.73000 

-.04207 

$ 

.74000 

-.04080 

$ 

.75000 

-.03950 

$ 

.76000 

-.03819 

$ 

.77000 

-.03687 

$ 

.78000 

-.03552 

$ 

.79000 

-.03417 

$ 

.80000 

-.03279 

$ 

.81000 

-.03140 

$ 

.82000 

-.02999 

$ 

.83000 

-.02856 

$ 

.84000 

-.02712 

$ 

.85000 

-.02566 

$ 

.86000 

-.02418 

$ 

.87000 

-.02269 

$ 

.88000 

-.02117 

$ 

.89000 

-.01964 

$ 

.90000 

-.01810 

$ 

.91000 

-.01653 

$ 

.92000 

-.01495 

$ 

.93000 

- .01334 

$ 

.94000 

-.01172 

$ 

.95000 

-.01008 

$ 

.96000 

- .00842 

$ 

.97000 

-.00674 

$ 

.98000 

-.00504 

$ 

.99000 

-.00332 

$ 

1.00000 

-.00158 

$ 


END* $ 


####################### ################################################## 

Generate Blade Tables 

############################# ##### #### ################################### 

Now generate the blade database and output to a library 





ROTO NET Primer 


('-8 


called DEM01.WRK; 


EXECUTE IBS $ 

EXECUTE IBA $ 

EXECUTE IBL $ 

$ 

S 

UNLOAD /DEMOO 1/ IBS, IBA, IBL, EDGE $ 

$ 

$ 

ENDCS $ 
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Example Deck 2 


ANOPP JECH0= . TRUE . JL0G=. FALSE. NLPPM=55 $ 

STARTCS $ 

$ 

$ fixes; 

$ activate ABL-STL-TIN executes 

$ need flight conditions for tie-down test 

$ check LRP->LRN data flow for changed parameters 

$ how to do observers for RBN? 

$ remove EDGE(THICK) , IBA(ZERO) for level 10 

$ 

$ ######################################################################### 

$ 

$ File; demo02.inp 

$ 

$ Author; Devon Prichard 

$ 

$ Level; 03/02/09 

$ 

$ Purpose; this input deck computes the performance, tone noise and 
$ broadband noise of the MD-500E main rotor. 

$ 

$ Input Libraries; LIBRARY. WRK, DEM00 l.WRK 
$ Output Libraries; DEM002 . WRK 
$ 

$ Modules used; 

$ 

$ APM - atmospheric properties module 

$ LRP - rotor performance module 

$ LRN - rotor tone noise module 

$ RBN - rotor broadband noise module 

$ ABL - atmospheric boundary layer module 

$ STL - streamline tracing module 

$ TIN - turbulence ingestion noise module 

$ 

$ ############################################################### ########## 

$ Setup Parameters 
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$ ######################################################################### 

$ 

$ Load atmospheric data to get sound speed, etc. 

$ 


PARAM IUNITS = 7HENGLISH $ use english units 

$L0AD /LIBRARY/ ATM=ATM0S (TM0D=STRD) $ load atmosphere from LIBRARY. WRK 
EXECUTE APM $ 

SPURGE ATM $ 

$ 

$ Define parameters unique to this rotor, otherwise use defaults; 

$ 

$ 


PARAM 

A0 

PARAM 

BETA 

PARAM 

F 

PARAM 

I OUT 

PARAM 

IMPR0V 

PARAM 

IPRINT 

PARAM 

MACHRF 

PARAM 

NBLADE 

PARAM 

NC0EFF 

PARAM 

NHARM 

PARAM 

NTIME 

PARAM 

NU 

PARAM 

R 

PARAM 

RX 

PARAM 

LBLCAL 

PARAM 

TB2CAL 

PARAM 

TBLCAL 

PARAM 

TEBCAL 

PARAM 

TETCAL 

PARAM 

TVFCAL 

PARAM 

WB 

PARAM 

ZETA 

EVALUATE RS 


$ flight condit 


= 0.2 

= 0.0 0.0 0.0 
= 0.0344 
= 1 

= .TRUE. 

= 1 

= 0.60 
= 5 
= 3 
= 8 
= 32 

= 0.0001572 
= 13.205 
= 5.0 
= .FALSE. 

= .TRUE. 

= .FALSE. 

= .FALSE. 

= .TRUE. 

= .TRUE. 

= 0.57E-04 
= 0.0 0.0 0.0 
= RX * R 


$ initial collective set {recall this 
$ is relative to root} in radians 
$ initialize flapping to zero 
$ vehicle flat plate drag area 
$ use spherical option 

$ use blade properties from IBS-IBA-IBL 
$ print input data 
$ rotor tip Mach number 
$ number of blades per rotor 
$ number of flapping/lagging coefficients 
$ number of desired harmonics 
$ number of desired time steps 
$ kinematic viscosity <3 sea level std., ft~2/s 
$ blade length in feet 

$ nondimensional sound bubble radius {re R} 

$ assume no LBL-TE noise 
$ use Brooks TE bluntness noise model 
$ don't use Schlinker-Amiet TBL-TE model 
$ don't use Grosveld TE bluntness model 
$ use Brooks TBL-TE model 
$ include vortex formation noise 
$ total blade weight 
$ initialize lagging to zero 
$ dimensional sound bubble radius 


$ 


PARAM 

EVALUATE 

EVALUATE 

PARAM 

EVALUATE 


KTS =95.0 $ 

VA = 1.69 * KTS $ 

MZ = VA / CA $ 

WT = 0.0051 $ 

RINF = CA * R / NU $ 


aircraft speed in kts 
flight speed in fps 
forward flight Mach number 
weight coefficient 

scaling Reynold's number for broadband 


$ 

$ 

$ ######################################################################### 

$ Setup Members 
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$ ################################################################ ######### 


LOAD /DEM001/ 


$ load blade database members from demol. 


$ Use finer spanwise grid than for blade properties run. 

$ 

UPDATE NEWU=GRID S0URCE=* $ 

-ADDR 0LDM=* NEWM=XI1 F0RMAT=4H*RS$ $ 

0.200 0.300 0.400 0.500 0.600 

0.700 0.750 0.800 0.850 0.875 

0.900 0.925 0.950 0.975 1.000 $ 

-ADDR 0LDM=* NEWM=PSI F0RMAT=4H*RS$ $ 


0 

.00 

0 

.05 

0 

. 10 

0 . 

. 15 

0 

.20 

0 

.25 

0 

.30 

0 

.35 

0 

.40 

0 

.45 

0 . 

.50 

0 

.55 

0 

.60 

0 . 

.65 

0 . 

.70 

0 . 

.75 

0. 

.80 

0 , 

.85 

0 . 

.90 

0, 

.95 

1 . 

,00 










END* $ 


TABLE TIP (ANGLE) 
INT = 1 
IND1 = RS 
DEP = RS 


S0URCE=* 


END* $ 


UPDATE NEWU=SFIELD SDURCE=* $ 

-ADDR 0LDM=* NEWM=THETA F0RMAT=4H*RS$ $ 


0.0 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

110.0 

120.0 

130.0 

140.0 

150.0 

160.0 

170.0 

180.0 

$ 

-ADDR 

0LDM=* 

NEWM=PHI 

F0RMAT=4H*RS$ $ 


60.0 -45. 

0 -30.0 

-15.0 

0.0 15. 

0 30.0 

ADDR NEWM= 

: FREQ 0LDM=* F0RMAT=4H*RS$ $ 

12.5 

16.0 

20.0 

25.0 

31.5 

50.0 

63.0 

80.0 

100.0 

125.0 

160.0 

200.0 

250.0 

315.0 

400.0 

500.0 

630.0 

800.0 

1000.0 

1250.0 

1600.0 

2000.0 

2500.0 

3150.0 

4000.0 

5000.0 

6300.0 

8000.0 

10000.0 

$ 


END* $ 


######################## ################################################# 

Compute Airloads and Noise 

########################### #### ## ######################################## 
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UPLIST $ 

EXECUTE LRP 
EVALUATE BETA(l) = AOL 
EVALUATE BETA (2) = AIL 
EVALUATE BETA (3) = B1L 
PARAM IPRINT = 3 
EXECUTE LRN 
PARAM I OUT = 1 

EXECUTE RBN 

EXECUTE ABL 
EXECUTE STL 
EXECUTE TIN 

UNLOAD /DEM002/ SFIELD , LRN, 

$ 

$ 

ENDCS $ 


$ do rotor performance calculation. 
$ put flapping into LRN form 
$ 

$ 

$ print input and output data 
$ do tone noise. 

$ output SPLs only 
$ do broadband noise. 

$ 

$ 

$ 


RBN $ 
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Example Deck 3 


ANOPP JECH0= . TRUE . JL0G=. FALSE. NLPPM=60 $ 

STARTCS $ 

$ 

$ modify NOISE load for TIN members 
$ modify PRO execute " " " 

##################### #################################################### 

File; demo03.inp 

Author; Devon Prichard 
Level; 03/02/09 

Purpose; this input deck calculates the aircraft flight trajectory 

and propagates the source noise for the MD-500E helicopter. 

Input Libraries; LIBRARY. WRK, DEM002.WRK 
Output Libraries; none 

Modules used; 

ABS - atmospheric absorption module 
SFO - steady flyover module 
GEO - source/observer geometry module 
APM - atmospheric properties module 
PRO - broadband propagation module 
PRT - tone propagation module 
LEV - noise level integration module 
EFF - EPNL calculation module 

#################### ##################################################### 

Setup Parameters 

###################### ################################################### 
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PARAM IOUT =1 $ make SFO output FLI(PATH) 

PARAM I UNITS = 7HENGLISH $ use english units 

PARAM IPRINT =3 $ print input and output 

$ 

$ define trajectory stuff 

$ 

PARAM JF = 200 $ total number of time steps 

PARAM XA = -2000.0 $ starting location 

PARAM XF = 2000.0 $ ending location 

PARAM YA = 0.0 $ starting location 

PARAM ENGNAM = 3HXXX $ make consistent with LRN, RBN defaults 

PARAM DELF =41.0 $ narrowband bandwidth, Hz 

$ 

$ flight conditions 

$ 

PARAM ALPHA = -0.0981 $ shaft angle from LRP calculation {rad} 

PARAM KTS =95.0 $ aircraft speed in kts 

EVALUATE VA = 1.69 * KTS $ flight speed in fps 

PARAM ZA = 250.0 $ starting location 

$ 

$ ########################################################### ############## 

$ Setup Members 

$ ############################################### ########################## 

S 

LOAD /LIBRARY/ ATM=ATM0S (AAC=SAEABS TM0D=STRD) $ 

LOAD /DEM002/ $ load source noise data from demo2 run 

$ 

UPDATE NEWU=0BSERV S0URCE=* $ 

-ADDR 0LDM=* NEWM=C00RD F0RMAT=4H3RS$ $ 


0. 

.0 

-500. 

.0 

0. 

.0 

$ 

0. 

.0 

-250 

.0 

0, 

.0 

$ 

0. 

.0 

0. 

.0 

0 

.0 

$ 

0 

.0 

250 

.0 

0 

.0 

$ 

0 

.0 

500. 

.0 

0 

.0 

$ 


END* $ 

$ 

S 

$ ######################################################### ################ 

$ Compute Trajectory Info 

$ ######################################################################### 

$ 

EXECUTE ABS $ atmospheric absorption 

EXECUTE SFO $ steady flyover trajectory 

EXECUTE GEO $ source/observer geometry 

$ 

$ 

$ Load atmospheric data to get sound speed, etc. 

$ 
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EXECUTE APM 
PURGE ATM 


Define parameters unique to this rotor, otherwise use defaults; 


PARAM I OUT = 1 
PARAM IPNLT = .TRUE. 
PARAM I PR I NT = 0 
PARAM MEMSUM = 5HPRD 
PARAM MEMSUMN = 5HPRT 
PARAM PROSUM = 4HRBN 
PARAM PROTIME = 3HXXX 
PARAM PRTTIME = 3HXXX 
PARAM SURFACE = 4HHARD 
PARAM R = 13.205 

PARAM RX = 5.0 

EVALUATE RS = RX * R 


$ output SPLs 

$ ask for PNLT calculations 
$ print nothing 

5HPRES $ have LEV use PRO results 
5HPRES $ have LEV use PRT results 
$ have PRO use the RBN tables 
$ have PRO use the RBN tables 
$ have PRT use the LRN tables 
$ use hard surface model 
$ blade length in feet 

$ nondimensional sound bubble radius {re R} 
$ dimensional sound bubble radius 


######################################################################### 

Propagate Noise 

######################################################################### 

EXECUTE PRO YYYYYY=RBN GEOM=BODY $ propagate broadband noise 
EXECUTE PRT YYYYYY=LRN GEOM=BODY $ propagate tone noise 
EXECUTE LEV $ compute OASPL, PNLT, etc. 

J ARAM IPRINT =2 $ print EPNL outputs 

EXECUTE EFF $ compute EPNLs 


ENDCS $ 
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APPENDIX F 


Glossary 


liequenoy weighting wit li higlicsl gain at around 
2000 Hz. Used lor approximating the varying 
liunian response to sound at differing frequencies 

loss of wave energy due to viscous losses, heat 
losses, molecular energy losses 

aerodynamic forces on the rotor blade, either 
expressed as surface pressure or lift per unit 
span 

a range of frequencies, typically used to 
partition the* total trequency range of interest 
into a set of adjacent bands 

Blade Passage Frequency, equal to 2; riLV/, 

Computational Fluid Dynamics, a category of methods 
used to compute the nonlinear aerodynamic response 
of complex shapes to the flow conditions 
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diffract ion 


I)NL 


Doppler shift 


emission time 


FPNL 


deflection of wave path around an obstacle 
which is large compared to the acoustic wavelength 

Day-Night Level, sometimes denoted La,i* 
expresses! in dH; an A-weighted time 4 average of noise 
oven* a 21 hour period, requiring knowledge of flights 
per hour. 

the change in frequency due to relative motion 
between the source and observer, for example 
the Doppler shift for a source moving direct ly 
toward a stationary observer is 

9 = (^) 

where; / is t he frequency of t he source when stat ionary. 
/' is the perceived source frequency with relative motion. 
\ s is the source 4 speed (positive 4 on approach) 

t ime 4 at which e 4 vent had to be emit ted by t lie' source 4 
in order for the 4 event t ej be 4 received at time 4 / 
at the observer location. Also called source* 
t ime 4 and ret aided time. 

(Infective Perceived Noise* Level) a noise metric 
based on time integration of PNLT 


FW-H e 4 e] nation 


the 4 Ffowcs Williams-llawkings equation. a 



HOTONKT Primer 


form of wave' equation governing aeroacoustic noise 
generated by an arbitrary body (see Equation d.l) 


he*]p fib's 


documentation for all Functional Modules 
in VAX/VMS Help format (also called "pseudocode") 


Hz 


Hertz, frequency in cycle's per second 


interference 


superposition of waves wit h same' or ne*ar frequencies 


ban see DNL 

log ide'iilitie's .v = a u => ;/ = Iog f( .r 

b>g( n/) = log .r + log// 
log ~ = log a* - log // 
log(.r /J ) = n log .r 

man page's documentation for all Functional Module's 

in Hnix man format 

narrowband band which is wider than 1 ( X of center frequency 

but smaller than ^-octave': for rotorcra.fi 
noise, bandwidth of narrowband analysis 
analysis ranges from 1 Hz to 10 Hz. 

OASPL (OverAll Sound Pressure' Level): a frequency 

integration of a noise' spectrum at a given i 
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observer location and reception time, 
with equal weighting for each frequency band 

observer time t ime at which a sound is heard or received, also 

called reception time 

Pascal denoted Pa. a unit of pressure in SI (N/m 2 ) 

PNL Perceived Noise Level 

a weighted frequency integration of noise* at 
a particular observer position and time*, 
similar to dBA in concept 

PNLT Perceived Noise' Level (Tone-Corrected) 

based on adding tone corrections to the PNL-weighte 

^-octave band spectrum, then 

performing a frequency integration as for PNL 

pitch pure tone frequency (standard acoustics usage). 

also the aerodynamic angle relative to the incident 
airstream (standard aerodynamics usage) 

reception time time at which a sound is heard or received, also 

c alled observer time 

reflection a return of wave energy partly in tlx* direction of 

the source, occurs at all interfaces in 
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refract ion 


retarded time 


seat 1 ('ring 


sound bubble 


Sound Pressure 1 


transmission media, (‘an be gas/gas, gas/liquid, 
gas/solid, etc. 

change in wave direct ion of travel due to change 
in medium; depending on the angle between wave' and 
interface, some energy passes through at a new angle, 
and some reflect. s. In aeroacoust ies. due to 
change in atmospheric' sound speed with altitude. 

time at which ('vent had to bo emitted by t he source 
in order for the ('vent to be received at time / 
at the observer location. Also called source 
t ime and emission time. 

deflection of wave path around obstacle' which 
is small compared to wavelength 

a hypothet ical sphere surrounding t Ik* aircraft 
where information is transform'd from the 
source noise models in the aircraft reference 
frame to the propagation models in the ground-fixed 
frame 

xnc! (SIM.) SPL = 20 log l0 (f t ) = 20 log l(J p + f) I (dH) 

= J01og 1u (jg ? ) = 10log 10 (// 2 ) + 91 (dH) 


source' t ime' 


time at which event had to be emitted bv the' source' 
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specific acoustic impedence 


speed of sound 


spherical spreading 


wavelength 


in order for the event to be received at time / 
at the observer location. Also called retarded 
time and emission time. 

SPL ]2 = 201o glo [clSPLt/m + r<.s •PL 2 im] 


Z = 

r 

( pressure/])art ide velocity ) 

( length /t ime) 

reduction in pressure wave strength with increasing 
time; for the spherical wave emanating from an ideal 
point source, the pressure field can be written as 
/>(/’./) = 

where k — u?/r. and A is a magnitude constant. 

AS PI. = 20 log )0 (pz/l’t ) = 20 log >) 

for r <2 = 2r\. ignoring harmonic component. 

=► ^SPL = 20\og w (\) = -(ulH 

i . e . a 6 dB reduction in SPL with ('very doubling of 
distance. 



distance between wave peaks, meters 



wave number 


number of waves per unit length. I /meter 
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wave equation 


yoclavr band 


r- i 


an equation describing phenomena with wave 
behavior, in this case noise, which has the form 



for the acoustic pressure field of an inviseid, 
ideal gas with constant sound speed 

interval between two frequencies having a \/2: 1 ratio, and 
successive (‘enter frequencies having a ratio 
of 5: | (standard center frequencies are 

50 Hz, 6:b 80. 100. ...) 
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Symbols 


(IS 

T 

f 

//( ) 

/ 

j 

h 

lj 

lr 

m 

\! 

Mj 

Mr 

A 

A'/, 

n 

P 

P 


sound speed, meters /second 

differential unit of body surfa.ee. in 2 

body function for which T — 0 defines the surface 

frequency. Hz 

Hoax i side step operator 

v^T 

coordinate index used for indicia! vector notation 
wave number. 1 /meters 

aerodynamic force in the /- tli direction. S/nr 
aerodynamic force* along tin* source/ observer direction. N / m 2 
direction index 
Mach number 

source velocity in the y - 1 h direction 
Mach number of source/observer relative 
motion in the radiation direction 
number of data points in a time series 
number of blades on tin* rotor 
direction index, band counter 
acoustic pressure. N/m 2 
total pressure, Pascals 
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P = P + P 

rol’crencp pressure, N/m 2 (2.0x10 ’ lor air) 
roof mean square pressure. N/m 2 


P 


PL 

PT 

(P 2 ) 

V 


I 

Tjj 

it 


r 


r 


a 


if 

A . V, Z 


A 


moan (time-averaged) pressure, Pa 

loading noise, N /m 2 

thickness noise. N/m 2 

mean square pressure. {N/m 2 } 2 

Fourier transform of acoustic pressure 

distance between source and observer, meters 

radiation vector, i . e . vector between source and 

observer, meters 

y-th component oi the radiation vector, 
normalized by r 

observer (reception) time, seconds 
Light-hill stress tensor 
particle displacement, meters 
particle velocity, meters/second 

local normal velocity ol blade surface, meters/second 

observer position vector, meters 

source position vector, meters 

ground-fixed coordinate's, meters 

(A’ is the longitudinal direction. 

) is the lateral direct ion. 

Z is the vertical direction) 
specific heat ratio 
wavelength, meters 
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p 

air density. kg/m" 

/'x 

freest, ream air density, kg/m'* 

T 

source (emission) time, seconds 

n 

rotor rotational speed, radians/second 

a 1 

circular frequency* Hz 


= '^f 

i) 

partial derivative operator 

H ) 

Dirac delta function 

V 

gradient operator 

( )„, 

quantity at retarded (('mission) time 

(') 

time derivative' 
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